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Abstract. This paper considers the question of global in time existence and asymptotic be- 
havior of small-data solutions of nonlinear dispersive equations with a real potential V. The 
main concern is treating nonlinearities whose degree is low enough as to preclude the simple 
use of classical energy methods and decay estimates. In their place, we present a systematic 
approach that adapts the space-time resonance method to the non-Euclidean setting using the 
spectral theory of the Schrodinger operator — A + V. We start by developing tools of indepen- 
dent interest, namely multilinear analysis (Coifman-Meyer type theorems) in the framework of 
the corresponding distorted Fourier transform. As a first application, this is then used to prove 
global existence and scattering for a quadratic Schrodinger equation. 
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1. Introduction 

1.1. The homogeneous setting. Consider a prototypical homogeneous nonlinear dispersive 
PDE of the form 

/ id t u + L{D)u = Af{u), D:=-iV , , 

\ u (t = 0) = u [LL) 

where the equation is set in R d , u is a complex-valued scalar or vector, L(D) is a real Fourier 
multiplier, and N is superlinear of order p in u (e.g., |it| p ). Taking L = A, for example, gives the 
nonlinear Schrodinger equation (NLS); L = \D\ corresponds to the wave equation; and L = (D) 
is the Klein-Gordon equation. We call (11.11) homogeneous due to the homogeneity of the linear 
part idt + L(D). We might also refer to it, variously, as the flat, unperturbed, or Euclidean 
problem. 

As far as the global theory of this equation goes, the most basic question to ask is: With uo 
sufficiently small, smooth and localized, does there exist a global solution? If so, does it scatter? 
For high degree nonlinearities, the answer is yes, almost regardless of the structures of L or J\f: 
dispersive or Strichartz estimates for the linear part L are enough to close the argument. In 
contrast, classical energy and decay estimates fail for low degree nonlinearities (p less than the 
Strauss exponent), and the structures of L and M, or more precisely resonances, start playing a 
decisive role. In fact, global existence is not assured, even for small smooth data (cf., e.g., |28j). 

As it turns out, many physically interesting problems (in general relativity, plasma physics, 
water waves, etc.) fall into the second category. Consequently, this regime has been studied 
intensively over the past forty years. Most notably, two main methods were devised in the 
1980s to deal with resonances: the normal forms approach of Shatah |40j . and the vector fields 
approach of Klainerman |31| . Recently, the first author, Masmoudi and Shatah introduced a 
new method based on the concept of space-time resonances [141 \T7\ [22] which brings together 
(and goes beyond in some cases) the normal forms and vector fields methods. 

1.2. State of the art in the inhomogeneous setting. In this paper, we seek to give a sys- 
tematic treatment to perturbations of (II. lj) by a time-independent potential, namely equations 
of the form 

/ id t u + L(V-A + V)u = N{u) . , 

\ u(t = 0) =u . { } 

Because of the appearance of the operator \J — A + V, we refer to this type of problem as 
inhomogeneous, distorted, perturbed, or non- Euclidean. 

There are many reasons, both mathematical and physical, to study equations of this form. 
First, in a general sense, it is important to understand how properties of solutions to the ho- 
mogeneous problem react to various perturbations, such as external forcing which is commonly 
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represented by introducing a potential. Second, even working with a homogeneous model, estab- 
lishing the stability of a bound state or traveling wave often requires studying an inhomogeneous 
problem. This is a simple consequence of the fact that the linearized operator around a ground 
state automatically involves a potential. Third, many important physical systems come with 
a potential inhomogeneity intrinsic to the model. For example, one encounters these issues in 
studying the stability of the equilibrium state for the water waves equation, since the resulting 
linearized problem will be homogeneous only if the bottom is perfectly flat [17]. In general 
relativity, as well, one often encounters nonlinear wave equations with a potential (e.g., when 
studying Schwarzschild metric in Regge- Wheeler coordinates [8].) 

Regarding (jl.2p . again the fundamental questions to ask are: For uq sufficiently smooth, 
localized, and small, does there exist a global solution? If so, what is its asymptotic behavior? 
As in the homogeneous setting, if the nonlinearity is of high enough degree, global existence and 
scattering follow easily. The only additional proviso is that the operator idt + L(\/— A + V) must 
enjoy similar dispersive estimates as the unperturbed operator idt + L{D). Sufficient conditions 
for this have been studied by many authors (cf., e.g., [38], [39] ) . 

On the other hand, comparatively little is known for low degree nonlinearities. A first line 
of research in this direction was pursued by Keel, Smith, and Sogge [30], Metcalfe, Nakamura, 
and Sogge [35], Metcalfe, and Sogge [36]. These authors considered the nonlinear wave equation 
in dimension 3 on an exterior domain, and assuming that the nonlinearity satisfies the null 
condition. Though this is not explicitly in the framework of (jl.2p . note that one can think 
of exterior domains as a limit of potential perturbations. In these works, small-data global 
existence was proved by following either of two distinct approaches: the conformal method 
of Christodoulou, or the vector fields method of Klainerman. To adapt the latter approach, 
the authors had to (smoothly) cutoff the standard full-space vector fields on a compact set 
containing the forbidden domain and prove weighted space-time energy estimates. Very recently, 
S. Yang [51] proved global existence for a nonlinear wave equation involving linear perturbation 
terms like L^d^u. This was done under either smallness or mild (time-)decay assumptions on 
by using the vector field method and a new scheme developed by Dafermos and Rodnianski [11] . 
We remark that none of these authors address the question of asymptotic behavior of solutions, 
which is something that we are able to do here. 

Turning to the Schrodinger equation, Cuccagna, Georgiev, and Visciglia [10] considered the 
nonlinear Schrodinger equation in dimension 1, with an added potential V. They were able to 
prove global existence and describe the asymptotic behavior by adapting the (one-dimensional) 
approach in [34]. All the works which we have discussed so far examine (as we will be doing 
for the rest of this manuscript) the situation where the linear operator spectrum is continuous, 
but one cannot get global existence simply by means of the dispersive estimates, since the decay 
rate given by the linear group is too weak and the power of the nonlinearity is not high enough 
to compensate for it. An interesting related problem is the situation where the strength of the 
dispersion is sufficient, but the potential exhibits linear, or nonlinear, bounds states. We refer 
for instance to [421 H5| I21T d] for more on this direction of research. 

1.3. Space-time resonances in the inhomogeneous setting. The central objective of this 
work is to develop a space-time resonance method in the inhomogeneous regime. At a conceptual 
level, this would yield a more general understanding of resonant structures in the distorted 
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setting. In a more practical sense, this will allow us to address both the global existence and 
the asymptotic behavior for a number of problems of the form (|1.2p . 
To make matters more concrete, let us consider the equation: 

{idtu — Au + Vu = u 2 
U[t = (J) = Uq. 

set on M 3 . Note that a quadratic nonlinearity qualifies as low degree in M 3 as the Strauss 
exponent there is exactly 2. In particular, this indicates that resonances need to be taken into 
account. We shall argue here in an informal fashion. We start by introducing the generalized 
eigenvectors e(x,£) which diagonalize the self-adjoint operator H = — A + V. They are given by 

He(x,£) = |£| 2 e(x,£) and e(x, £) ~ e ix< if |x| ->• oo. 

The properties and necessary conditions for the existence of these generalized eigenfunctions are 
reviewed in Section 12.11 They naturally give rise to a distorted Fourier transform, which we 
denote by J 7 " or 



(J*9)(t) = <?*(£) := T^L^ / 9(x)e(x,0dx. 



(2vr) 3 / 2 

(cf. Theorem 12. ip . As in the homogeneous setting, the idea is to write Duhamel's formula for 
the profile / := e~ ltH u in the distorted Fourier space. This gives after a short computation 



/»(*,£) = «S(£)-i f JJ e- ls ^^'fi(s,v)fKs,C)M(^ V ,C)dvdCd S (1.4) 
where 

M(Z,7 1 X) = j^^ J e(x,Oe(x, V )e(x,()dx and 77, C) = |£| 2 + |r/| 2 + |C| 2 - 

Observe that the flat case, M. (£, 77, C) = <5r3 (^ + 77 + £), and thus the above integral reads 



This simpler formula is of course a manifestation of the convolution identity fg = f*g, which 
does not hold anymore for the distorted Fourier transform. As we argue below, this is the source 
of tremendous difficulties in the analysis if V 7^ as it makes manipulations on the (distorted) 
Fourier side much more delicate. But it is also responsible for an interesting new phenomenon: 
all frequencies £, 77 and £ interact, not only the ones that add up to zero. However, we expect 
Ai to be nicer away from the hyperplane {£ + 7/ + £ = 0}; on this set, a singularity occurs, which 
seems to be in general more complicated than the <5-function we see in the flat case. 

We now give a brief outline of the space-time resonance method (see |14| for a more complete 
discusssion in the case V = 0). The idea is to consider the integral on the right-hand side 
of (II. 4j) as an oscillatory integral. Proving global existence and scattering essentially amounts to 
ensuring that this integral remains bounded (in a sense that we keep vague for the moment) as t 
goes to infinity. In the flat setting, obstructions to this behavior come only from the stationary 
points (in the (s, 77) integral) of the phase <j> in (|1.5j) . This leads us to define the time and space 
resonant sets as the sets where <f> and d v (p vanish, respectively. In contrast, in the distorted 
setting (|1.4p . there are three integration variables (s,7],() and a distribution M. that is singular 
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on the set {£ + f? + C = 0}> a fact which poses a challenge even in defining the space resonance set 
in the first place. More importantly, it is crucial to point out that there are specific directions 
(or vector fields) along which one can differentiate M without increasing the severity of its 
singularity (as quantified by its order as a distribution, sajQ). This phenomenon is even lurking 
in the flat case: applying the vector fields 8^ — 8^ to <5(£ + rj + £) gives a distribution of order 
0, whereas applying d^a yields a distribution of order 1. 

Understanding and coming to terms with these "admissible" directions plays a central role 
both in defining the space resonant set below, and, more generally, in bounding expressions 
involving derivatives of the distribution M. In Section 13.41 we show that, similar to the case 
V = 0, the directions given by the vector fields 8^* — d^a are, morally speaking, admissible. 
More precisely, though, the analysis suggests the need for certain non-local modifications of 
those vector fields (cf. Section [3.4.1 1 and identity (|3.56p ). but let us gloss over this important 
fact for now. 

With this information in hand, one can finally mimic the flat scenario and make the following 
definitions. 

• The time-resonant set T consists of points that are stationary in s (i.e. 8 s (s(ft) = 0): it 
is thus T = {(ft = 0}. This corresponds to resonant interactions in the ODE sense. 

• The space-resonant setS consists of points that are stationary in {n, Q when we differen- 
tiate the phase function along the admissible directions. Therefore, S = {{d^—d^cft = 0}. 
This corresponds to wave packets traveling at the same group velocity and agrees with 
the previous definition when V = 0. 

• Finally, the space-time resonant set 1Z is the intersection of the space and time resonant 
sets: K = T D S. 

This classification of stationary points gives a clear intuitive understanding of the resonant 
interactions at hand, and this understanding helps to answer the next question: How can we 
obtain the required estimates? 

• Away from the time resonant set, it is possible to integrate by parts in s using the 
identity j^8 s e ls ^ = e ts ^ . This corresponds to a normal form transform, which effectively 
increases the degree of the nonlinearity. 

• Away from the space resonant set, it is possible to integrate by parts in (rj, C) using 
the identity i8(9 llg ^ (fy" - 8 c <x)e is< t> = e is ^ . Of course, the gain here comes from the 
factor of s in the denominator which improves the convergence chances of the integral. 
Note that the use of an admissible vector field is essential to control the terms involving 
derivatives of A4 which are inadvertently produced when we integrate by parts. For the 
same reason, we require bounds on derivatives of or, equivalently, weighted-norm 
estimates on /. 

• Finally, there remains the space-time resonant set, for which none of the above strategies 
applies. Here the analysis becomes more problem-dependent. One possibility would be 
to simply split the above integral by cutting off a shrinking neighborhood of the space- 
time resonant set. To estimate the piece thus removed, one can use its shrinking size to 
gain integr ability; as for the rest, the above integrations by parts apply. 



Recall that a distribution F is said to be of order m if (r, <fi) can be bounded in terms ||</>||c m for any test 
function <j>. Distributions of order 0, like M, are measures |24| . 
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For equation (II. 3p . a simple computation shows that 1Z = T = {(0,0,0)}, and thus space 
resonances can be ignored for this problem. Nonetheless, even in such a simplified context, a 
quantitative understanding of the "admissible derivatives" of the distribution A4 is indispensable 
if one needs to prove weighted estimates on the solutions. We will treat equations with S (jLT 
in a future work using the harmonic analytic framework developed here. 

1.4. Assumptions on the potential V . Before stating our main results, let us now describe 
the assumptions we shall impose on V. 

1.4.1. A qualitatative description. 

• HI Existence of distorted Fourier Analysis [1] 

• H2 Absence of discrete spectrum for — A + V 

• H3 LP boundedness of the wave operator fi := lim^-oo e itH e~ itA [SI 147]. 

In total, HI, H2, H3 amount to some regularity and decay requirements on V. For fur- 
ther discussion of HI and H2, including explicit sufficient conditions, see Section 12.11 (namely 
Theorem 12.11 and Remark |2.2jl . Assumption H3 has been proved to hold for 1 < p < oo in 
all dimensions d > 2 for so-called generic-type potentials V (i.e. with no null resonances) with 
sufficient decay at infinity. The literature on this topic is quite lengthy, and we give a somewhat 
detailed account in Section [2.21 

At various points in our analysis, we will need to take a deeper look at the wave operator £1. 
In those cases, we will assume one or both of the following: 

• H3* W 1,p boundedness of the wave operator Q [6j [47] 

• H3** An explicit structure theory for fi [6| I47j 

Each of these is closely related to H3; in particular most proofs of H3 and H3* rely precisely 
on an explicit structure theory as assumed in H3**. The reader is again directed to Section \2. 21 
for a more detailed presentation. 

1.4.2. Quantitative description on R 3 . On M 3 (where we set our PDE application), Beceanu 
provides in [6j the most refined result guaranteeing assumptions H3, H3*, H3** with minimal 
decay conditions on V (e.g., V £ (x)~ l ' 2 ~ e L 2 for H3, H3**). However, at one juncture in our 
analysis (Section 13.3ft we find it easier to use Yajima's slightly more explicit structure theory for 
Q in [47] , This comes at the cost of imposing more decay conditions on the potential in Theorem 
11.21 below (cf. Remark 12 .71) . A very simple condition (but very far from optimal) which ensures 
that all these hypotheses are met on R 3 is the following 

• The operator — A + V does not have discrete spectrum or null-resonances. 

• There exists e > such that 

\V(x)\ < {x)-^-\ 

1.5. Obtained results: Multilinear distorted Fourier analysis. Our first set of theorems 
lays the harmonic analytic groundwork for the space-time resonance method. We hasten to point 
out, though, that they are of interest in their own right. While it may seem that transitioning 
from the flat to the distorted regime is straightforward — all one must do is exchange J- for 

- in practice, it requires some new ideas, and a few highly technical arguments. In that sense, 
these theorems should be thought of as one of the principal contributions of this paper. 
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The oscillatory integrals encountered in Section 11.31 lead one to study multilinear operators 
of the type 

TU^'^jf 1 [[ m((,V,0fHv)9HC)M&vX)d(dri. (1.6) 

JJR d x'R d 

These are the generalizations of the classical pseudo-product operators of Coifman and Meyer [9] 
to the inhomogeneous setting. Typically, one would like to prove Holder-type estimates, i.e. that 
T is bounded from L q x LP — > LP' with p = — h ~ which hold trivially if m = 1. A natural 
condition on m mimics the Coifman-Meyer class of symbols by requiring 

|3£ • • • fl£m(fc, ...,*„)|< C ai ,..., a Mi\ + ■■■ + |£ re |)-(W+-+W). (1.7) 

Any m satisfying the above homogeneous bounds is said to be a Coifman-Meyer symbol (in n 
variables). 

Our first result is the distorted analog of the celebrated Coifman-Meyer theorem [9]. 

Theorem 1.1. Let V £ L d / 2 (R d ) (d > 2) be a potential satisfying HI, H2, and H3. Sup- 
pose that m = rn(^,rj,() is either a Coifman-Meyer symbol in three variables as in (|3.ip . or 
m (?> r ?>C) = m o( ? ?)C) where m,Q is a Coifman-Meyer symbol in two variables. 

(i) For any p,q,r £ (1, oo) satisfying - + - + £ = 1, we have 

Wr(fi9)\\ L r>( R d) < m ,V \\f\\Li(R d )\\9\\LP(R d ) (I- 8 ) 

provided that the Riesz transform = V(— A + V)" 1 / 2 is bounded on L p ,L q , and LP' . 

(ii) Suppose instead that V satisfies assumption H3*. Lf q,p,r,q,p G (l,oo) satisfy p = 

p + | = | + | _e / or some e > °' 

\\ T (f>g)\\L r '(M. d ) ^m,V,e \\ f \\ Li(R d ) llff IIlp (M d ) + II/IIl9(R' 1 ) Ilffllz/P (M d ) ) (I- 9 ) 

whenever f £ L q f) L q and g £ LP ' P\ LP '. 

The boundedness of the Riesz transform associated to a potential V on L p (M. d ) is clearly 
equivalent to the embedding of the distorted Sobolev space 

a question that has been an active subject of research for some time. We touch on the issue 
briefly in Section f2.4l However, for our PDE applications, the imperfect estimate (|1.9p suffices, 
and so we do not need to make such strong assumptions on V. The proof of Theorem 11.11 
starts by following the same strategy as that of the classical Coifman-Meyer theorem. This 
gives the needed bounds on the diagonal interactions (where the highest two frequencies are 
comparable) once one proves maximal and square function estimates for the distorted setting 
(cf. Lemma l3.3j) . In contrast to the flat setting, however, non-diagonal interactions are present 
and very problematic when there is a potential. To deal with them, we use a correlation identity 
involving three generalized eigenfunctions (the same idea was used in [23j for similar reasons). 
What results is a trilinear term involving the Riesz transform; this is precisely why we must 
require that £H is bounded on L p , or else sacrifice some e of integr ability. Related work, dealing 
with a more restricted class of operators, namely para-products instead of pseudo-products, can 
be found in [7]. 

Obtaining Holder-type estimates for T(/, g) turns out to be merely part of the story. Indeed, 
to successfully prosecute the space-time resonance method, we will need weighted estimates for 
T(f,g), or, equivalently, LP x LP — > LP estimates for (|1.6h with a derivative falling on the 
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integral. This leads inexorably to the issue of admissible directions discussed in Section 11.31 
Ultimately, we are able to obtain the needed estimates bounding \\(x)T(f,g)\\ Lr i in terms of 
IK^/lli' an d \\g\\LP f° r symbols that are slightly more regular than Coifman-Meyer ones and 
for potentials satisfying hypothesis H3**. This is the content of Theorem 13.131 which we do not 
transcribe here for brevity. 

1.6. Obtained results: Global existence for a quadratic Schrodinger equation. With 
the above tools in hand, we are able to prove the following. 

Theorem 1.2. Let V be a potential satisfying assumptions HI, H2, H3*, and H3** and 

consider the quadratic nonlinear Schrodinger equation in IR 3 (II. 3|) . Let X be the space defined 
by the norm 

WfWx := WfW^Hi + Wd^W^q- 

There exists eo > such that, for any initial data uq with \\u$\\ H i + Wd^u^W^ < eo, there exists 
a solution u £ X to (jl.3p defined for all time t > 0. In particular, 

\\<t)\\ L v < — forall2<p<6. 

t* ( p> 

Finally, u scatters in L 2 , namely e lt ^u(t) has a limit in L 2 as t — > oo. 

The proof of this result is presented in Section H] and is meant as a first application of the 
"inhomogeneous" space-time resonance strategy outlined in Section 11.31 The energy space X 
is natural in light of the dispersive estimates for the linear propagator e ltH (cf. Proposition 
I2.9p . In particular, we note that in the flat setting, global existence and scattering was recently 
proved in a similar space by Laillet [33], who simplified the argument based on the space-time 
resonance approach originally given in [16j . 

Finally, we remark that our choice of NLS with a u 2 nonlinearity is meant to provide a 
somewhat generic example of how to proceed when the unperturbed homogeneous problem can 
be treated using the normal forms approach. Many other homogeneous equations fall into this 
category, for example Klein-Gordon [lUJ [32] , water waves [SI H7] , Klein-Gordon-Zakharov [37] , 
and Euler-Poisson for the electrons |19} [T5] or the ions [20] . To see the full strength of the space- 
time resonance method, though, one must consider problems where the normal forms approach 
on its own is not enough. This we do in a forthcoming paper. 

Plan of the article. In Section [2j we present several known results on the spectral theory 
of Schrodinger operators —A + V which are needed in the rest of the paper. In Section 3, 
we develop the theory of multilinear distorted Fourier multipliers, that is to say the analog of 
pseudo-products when the distorted Fourier transform is substituted to the Fourier transform. 
This is applied in Section 4 to prove global existence and scattering for a quadratic Schrodinger 
equation in K 3 . 

Notation. We write A < B to signify that there is a constant C > such that A < CB. We 
also write A ~ B when A < B < A. If the constant C involved has some explicit dependency, 
we emphasize it by a subscript. Thus A < u B means that A < C(u)B for some constant C(u) 
depending on u. In some instances, we use the notation A ^ B to signify that the implicit 
constant C is large. Also, we use the Japanese bracket convention where (x) := \J\ + |x| 2 . 
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2. Tools from linear distorted Fourier Analysis 

In this section, we review a number of important topics from the theory of distorted Fourier 
analysis that will form the basis for our work in the remainder of the paper. 

2.1. Schrodinger operators. For a given potential V : R d — > R, consider the associated 
Schrodinger operator H := —A + V. Of particular interest is the situation where H is a 
perturbation of Hq := —A. For instance, when V £ L 2 (R d ), H can be realized as a self- 
adjoint operator on L 2 (R d ) with domain D(H) = H 2 {R d ). We may then ask: For what V do 
the spectral properties of H resemble those of HqI A natural starting point is to impose a 
compactness condition on the multiplication operator associated with V . With that in mind, 
we say that V is short-range (or, of class SR) provided that 

u G H 2 (R d ) H- (1 + \x\) l+e Vu E L 2 x (R d ) is a compact operator, 

for some e > 0. It was shown by Agmon (cf. pQ) that, for V of class SR, cr(iJ) = {Aj} Jg j U 
[0, oo); the continuous spectrum being [0, oo), and the discrete spectrum consisting of a countable 
set of real eigenvalues {Xj}, each of finite multiplicity. Furthermore, we have the orthogonal 
decomposition 

L 2 {R d ) = L 2 ac (R d )® Ll{R d ), 

where L 2 (R d ) is the span of the eigenfunctions corresponding to the eigenvalues {Xj}, and L 2 &c 
is the absolutely continuous subspace for H. 

It should be noted that numerous versions of this result predate the work of Agmon, only 
for different classes of potential. Indeed, pQ is a particularly significant milestone in a long and 
continuing effort to discern the optimal conditions on V under which these and other spectral 
properties of H hold (cf. the references contained in [l] for a summary of earlier works, and, 
e.g., [26] for an important recent improvement). 

In order to identify and study the resonances in PDEs like (jl.lj) or (]1.3p . it is best formulate 
them in frequency space. For this, we need a "well-behaved" eigenfunction expansion that 
diagonalizes H (on L 2 C ). Of course, this is in turn predicated on the existence of a natural 
class of generalized eigenfunctions that serve as analogs of the plane waves e lx ^ for Hq := — A. 
For each £ £ R d \ {0}, we know that |£| 2 is in the continuous spectrum of H; the associated 
eigenfunction is the distorted plane wave e(-;£) defined as the solution of 

ffe(.;0 = KI 2 e(-;fl, (2.1) 

with the asymptotic condition 

v(x;£) := e(x;£) - e lx '^ = 0(|x| _1 ) as \x\ ->■ oo, 
and the Sommerfeld radiation condition 

lim r(d r — i\£\)v = 0. 

r— >oc 

This can be expressed in a more convenient way via the resolvent: for z G C \ o~(H), define 
Ry(z) := (H — z)~ l and consider 

R7/(z) := lim Ryiz — ie). 

V e->0+ 
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The limit here is taken in the uniform operator norm topology on C((x)~ s L 2 , (x) s H 2 ) for s > 
2~, its convergence is the so-called limiting absorption principle (cf. [U Theorem 4.2]). The 
eigenfunction problem (|2.ip can then be recast as the Lippman-Schwinger equation: 

e(-;0 = e c - Ry(\^\ 2 )Ve^ e^x) := e ix <. (2.2) 

It can be shown that there exists a unique solution to (12.2() for any £ G IR d \ {0} provided that 
V = 0(|x| _1_e ) as |x| — > oo, for some e > (cf. |T|). Under this assumption, the distorted plane 
waves are relatively smooth in x, but have very little regularity in £. More precisely, for fixed 
(GK rf \{0}, 

e(-;0 G (z) s #|, for any s > (d + l)/2, (2.3) 

however, the map (x,£) \— > e(x;£) is merely measurable. One can improve this by requiring 
additional decay and regularity of V (cf., e.g., [25]). 

In view of the Fourier transform, we expect that the family {e(-;£)} forms a basis for the 
absolutely continuous subspace of H. This is indeed true, as was first proved by Ikebe [25] and 
later generalized by several authors. For consistency of presentation, we give here the version 
due to Agmon (cf. [TJ Theorem 6.2]). Before that, let us now impose assumption H2, namely 
that H has no discrete spectrum. However, we remark that many results in this paper (especially 
those in Section [3]) can be directly generalized to potentials with discrete eigenvalues by simply 
projecting on the absolutely continuous subspace L 2 C throughout. 

That said, the result is the following. 

Theorem 2.1 (Ikebe, Alsholm-Schmidt, Agmon). Consider the Schrddinger operator H with 
potential V satisfying H2 and 

(1 + |^|)2(i-t- e ) [ \V(y)\ 2 \y - x\~ d+d dy G L™(R d ), for some e > 0, < 9 < 4. (2.4) 
Jb 1 (x) 

Define the distorted Fourier transform F^ by 

(Pf)(0 := fKO ■■= jAjH P lim [ <^0mdx, (2.5) 

(2TT) d / 2 R-+ooJ Br 

where Br is the ball or radius R centered at the origin in R d . ThenF^ is an isometric isomor- 
phism on L 2 (W d ) with inverse formula 

f(x) = (p- 1 ft)(x):=-\- lim f e(x;0f*(Qdt. (2.6) 

Moreover, F^ diagonalizes H in the sense that, for all f £ H 2 (M. d ), 

Hf = F^MF^f, (2.7) 
where M is the multiplication operator u t— > \x\ 2 u. 

Remark 2.2 (Sufficient conditions for assumptions HI and H2). We are now, at last, able 
to give a precise meaning to assumption HI: we say that HI is satisfied provided that (i) the 
family of eigenf unctions {e(-,£)} exists with the regularity stated in (|2.3p . and (ii) the operator 
J 7 " defined by (|2.5p exists and exhibits the properties described in Theorem \2.1\ It follows that 
sufficient conditions for HI are that V satisfies H2, (|2.4|) . and V = 0(|x| _1 ~ £ ) as \x\ — > oo, 
for some e > 0. If we remove assumption H2, of course, we require only that F^ be a unitary 
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partial isometry with range L 2 c (M. d ). Note that, by imposing (|2.4p . we rule out the existence of 
nonnegative eigenvalues. In order for H2 to hold, we must require additionally that there are 
no negative eigenvalues, which is guaranteed, e.g., if the negative part ofV is not very large (for 
example if d > 3, Hardy's inequality implies that the condition V > — (d — 2) 2 /4\x\ 2 is sufficient 
to rule out both non-positive eigenvalues and resonances at as defined in (I2.1ip below). 

2.2. The wave operator 0. One consequence of Theorem 12.11 is that H and Hq are unitarily 
equivalent. To see this, note that by (|2.7D . 

H = nH n*, with Q := P~ X JF. 
The operator is called the wave operator. It can alternatively be defined by 

Q := s-lim e itH e- itH ° (2.8) 

t— >— oo 

in the strong operator topology. Note that some authors denote this as fi+, with J7_ being the 
result of taking the limit t — > +oo. Under the assumption that V is of class SR and H2 holds, 
these limits exist, and SI is a unitary operator on L? . This fact is often referred to as asymptotic 
completeness; for potentials satisfying (|2.4p . it is originally due to Agmon (cf. [TJ Theorem 7.1]). 
For us, the importance of ft lies in the intertwining relations 

e itH = n e itHo n*, Pn = F. (2.9) 

In other words, allows us to translate back and forth between the flat and distorted cases. 
Clearly, then, information about the structure and boundedness properties of £1 is extremely 
valuable. The foundational work in this direction is due to Yajima, who first proved the W k,p 
boundedness of f2 and Q* , under the assumption of sufficient smoothness and decay for the 
potential. We paraphrase his results below. 

Theorem 2.3 (Yajima, |47 t l4"9 l [T3" l I50|). Let k E N and consider the Schrddinger operator H 
with real potential V : M. d — > R for d > 3. Fix po, &o as follows: 

( po = 2, ko = if d = 3 

\ Po >d/2, k := L(d-1)/2J i/d>4. 

Assume that for some 5 > (3d/ 2) + 1, 

(x) <5 ||a Q T/|| L Po ( | a; _ y |< 1) G L™(R d ), for all a with \a\ < k + ko. (2.10) 

Then V is of class SR and so and Q* are well defined as operators on L 2 (]R rf ) n W k,p (M. d ). If 
we additionally assume that V is of 

Generic-type: there is no u G (x) e L 2 x (JL d ) solving Hu = 0, for any 6 > 1/2. (2-H) 

Then Q, and Q* may be extended to bounded operators defined on W k,p (M. d ). 

Remark 2.4. The assumptions here can be weakened somewhat when V is small or nonnegative. 
For example, in any dimension d > 3, hypothesis (|2,10p can be replaced by the following: for all 
\a\ < k, 

\\J r (x) a D a V\\ d -i is sufficiently small for some a > 2(d -2)/ (d - 1). (2.12) 
See \49\ Theorem 1.1, Theorem 1.2] and the surrounding remarks for further discussion. 
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Condition p. IIP asks for the absence of resonance at zero. It is indeed necessary for W k,p 
boundedness, though other forms of boundedness can be salvaged without it (cf., e.g., [48} I50j). 
We should also note that a similar result also holds on M 2 [27j . 

Recently, Beceanu [6] was able to extend a number of Yajima's results on R 3 , reaching scale 
invariant class of potentials. Namely, he was able to prove boundedness of the wave operator 
on L P (M. 3 ) for a class of potentials B such that (x)~ l l 2 ~ e L 2 C B C L 3 / 2 ' 1 . He also proved the 
W ,P (N?) boundedness of Q under the assumption that: 

( B p<3 

V £ < BnI 3 / 2+f p = 3 (2.13) 
[BnLP p>3. 

See [6] Corollary 1.5]. 

The above mentioned proofs are based on an asymptotic expansion of the wave operator as 
a Born series (obtained by repeated application of Duhamel's formula). The main difference 
between Yajima's and Beceanu's approaches is in bounding high order terms in this expansion. 
While Yajima resorts to a direct computation to estimate their contribution (see Lemma l2.5p . 
Beceanu uses an abstract version of Wiener's theorem which allows him to work under much 
lower decay assumptions on the potential V. 

Both Yajima and Beceanu prove their boundedness results by obtaining an "explicit descrip- 
tion" of the wave operator. We will need such a description in Section 13. 3} where we study 
commutators of £1 and position operators in IR 3 . We elect to use Yajima's version, which is 
slightly more explicit than Beceanu's. In effect, this requires us to impose Yajima's more re- 
strictive conditions on V in Theorem 13.41 (and hence Theorem 13.13} and Theorem 1 1 . 2[> . We 
believe that the weaker assumptions of Beceanu are sufficient, but showing this would entail 
considerable additional technical work. 



Lemma 2.5 (Structure of the Wave operator [58, 6]). Under the assumptions of Theorem \2. 
the adjoint ft* of can be written as: 

U*f = f-W 1 f + W 2 f - W 3 f + Lf (2.14) 

where W\, W2, W3 and L are bounded operators on L p for all 1 < p < 00 and have the following 
form: 

Wjf = / Kj(ti, ■ ■ ■ , tj-i, t, a/i, . . . , cjj)f(x + p)dti . . . dtj-idrdwi . . . dojj (2.15) 

J[0,oo)J- 1 xIx& 

Lf= f L(x,y)f(y)dy. (2.16) 



Here S = S 2 denotes the unit sphere in ]R 3 ; / = (—00, —2ujj(x + t\uii + . . . + £j-_io/j-_i)) is the 
range of integration of the variable r; and we denote furthermore 

y :=y- 2(u j .y)uj j , p := t\U)i + ...+ tj^iuJJZi - TUJj. 
The kernels K\, K2, K3, and L(x, y) are described more explicitly in Section [3.3\ and they satisfy 



l- K illi 1 ([o,oo)^,L 1 (s n )) ~v 1, 



sup / \L(x,y)\dx+ sup / \L(x, y)\ dy < v 1. 
yeR 3 JR 3 zeR 3 "/R 3 



(2.17) 
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Remark 2.6. It is easy to see from the above explicit description (particularly estimate (j2. 17j) ) 
that the operator Q is point-wise majored by a positive operato^ that is bounded on L P (M. 3 ) for 
all 1 < p < oo . This fact is also ( even more ) evident in Beceanu 's explicit description of £1 in 
[6], so we shall assume it throughout our work. 

Remark 2.7 (Sufficient condition for assumptions H3, H3*, and H3**). We say that H3 
holds provided fi is bounded on L p (M. d ), for any 1 < p < oo. Sufficient conditions for this are 
V € B (considered by Beceanu) and the generic-type assumption (|2.1ip when d = 3, or more 
generally (|2.10p and (|2.1ip with k = 0, for d > 4. (See also Remark \2.1$ . H3* holds provided 0, 
is bounded on W 1,p (M. d ). In M 3 , the optimal known conditions implying this are given by (|2.13p ; 
for other dimensions, the hypotheses of Theorem \2. 3\ with k = 1 seem to be the weakest currently 
available. We say that H3** is satisfied provided O can be written as a Born series of the form 
(|2.14p . This is satisfied whenever the hypotheses of Theorem \2.3\ holds (with k = 0). 

2.3. Distorted Fourier multipliers. Suppose V satisfies the hypotheses of Theorem 12.11 
Then, for any function m : M. d — > C, we define the distorted Fourier multiplier m{D«) to be 
the operator 

mp 1 ) -j'^m^J 1 , (2.18) 
where m(£) denotes the multiplication operator u \- > m(^)u. This is an analogue of the well- 
studied Fourier multipliers m(D) given by 

m(D) := F~ l m{£)F. 

The intertwining relation 

m(£> tt ) = ftm(V)fT (2.19) 

makes the boundedness of m(D^) on a given space equivalent to that of m(V) as soon as the wave 
operator is bounded. This immediately gives crucial boundedness properties for some standard 
(distorted) Fourier multipliers. 

For instance, we may define the (distorted) Littlewood-Paley operators as follows. Pick a 
smooth function ^ supported on -62(0), and equal to 1 on i?i(0). Let <I> := ^(2-) — ^. The 
Littlewood-Paley operators are then given by: 

Pn = ®{jP) and p<n = *(jP) Ne2 ^- 

Sometimes, we simply denote /at instead of -Pat/ so that the decomposition of a function / in 
dyadic frequency pieces reads 

/ = E p *f = E fa ( 2 - 2 °) 

Finally, if M is not a power of 2, it can be rounded to a power of 2, say N, and we set Pm = Pn- 
We also consider the distorted Sobolev spaces W^' p {R d ) := \D$\~ S L p x (R d ) (for -00 < s < d) 
and W^ p (R d ) = (Z? s ) _s Lg(M d ) (for all s G R). Thanks to the intertwining property and part 
(c) of Lemma 13.31 below, these spaces inherit the square function characterization from their 
Euclidean counterparts: for 1 < p < 00, 

w;- ~d ||(E ^ 2s I^WI 2 ) 1/2 IIl ? , \\f\\w r ~d||(E ^ 2s l^v/| 2 ) 1/2 || iS . 

JV62 Z 7VS2 N 



Recall that a positive operator is one such that Tf > whenever / > 0. 
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Here we adopt the convention that when Pi appears in a sum over dyadic blocks with N £ N, it 
should be understood as P<i- The question of the equivalence (or not) of the distorted Sobolev 
spaces and their homogeneous counterparts is addressed in Section 12.41 

The following proposition will be useful throughout the paper and its proof is direct by using 
the intertwining property (|2.19p . 

Proposition 2.8. Assume H3 is satisfied. 

(a) (Distorted Bernstein inequality) Let s > be given. Then, for each 1 < p < q < oo, and 
/€ 



lA/v^/llL9(K d ) + \\ P <l/y/if\\L9(K.*) £ s ^ 



LP 



(b) (Distorted Sobolev inequality) Ifl<p<q<oo and - — - = ^ , then 

iii^r/iu*<ii/iu* n/iiw < \\f\\wr- 

(c) (Distorted fractional integration) For each a > and t > 0, define the (distorted) 
multipliers 

A-<* :=t a / 2 <p(y/i\Dt\) a , 
where if is a smooth function satisfying 

l/k, k>2, 
1 k < I. 



<p(k) 



^/l<P<(7<oo and | — | < § orifl<p<q<oo and - — | = ^, i/ien 

l|A t - a /||w<^^ ( 5-^||/|| iP . 

Finally, we will sometimes need to use this inequality in conjunction with the linear 
group: under the same conditions as above, and assuming furthermore q < 2 < p, 

\\K a e~ aH f\\ Lq <t^T~\\\f\\ LP . 



2.4. Boundedness of the Riesz transforms. The Riesz transform associated with the elliptic 
operator — A + V is given by 

W/ := V(-A + VY^f =: ^Lf (2.21) 

The question of boundedness of the operator above has been the subject of intensive study in the 
past twenty years [HI [31 [H [18] as it is equivalent to boundedness of the embedding W$ ' P W 1,p . 
It is the analogue of the well-known question originally raised by Strichartz concerning the 
L P (M) boundedness of Riesz transforms , v on complete Riemannian manifolds (M,g) |44| . 

The above mentioned references give certain classes of potentials for which the Riesz transform 
(|2,2ip is bounded on various L p (R d ) spaces. We are not concerned here with the state of the art 
in that field, though, since we choose not to impose these conditions on V (especially that for 
a large class of potentials of interest to us, the operator 9\ is only bounded if 1 < p < d [18J). 
In fact, for the purposes of our PDE applications, we are content with the "imperfect" estimate 

clu. 
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Finally, note that if we regularize 9\ at low distorted frequencies, we directly get a bounded 
operator on L p (R d ) for all 1 < p < oo under our assumptions on the potential. Indeed, for V 
satisfying H2 and H3*, the operator 

33 : / V(/ - A + V)~ 1/2 f = V(^)" 1 / (2.22) 

is bounded on L p (E. d ) for all 1 < p < oo. This follows directly by noticing that (D^)^ 1 = 
il(V) _1 il*, using assumption H3* and the boundedness of V(V) _1 . 

2.5. Dispersive estimastes on M 3 . Dispersive estimates for the semigroup of a Schrodinger 
operator have been studied extensively; let us mention in particular the work of Journe, Sogge, 
and Soffer [29] , and the more recent papers of Rodnianski and Schlag [38] , and Beceanu [5] . A 
recent survey is given in |39j . 

Observe in any case that the L p boundedness of the wave operator Q corresponding to H 
implies that the group e ltH enjoys the same dispersive and Strichartz estimates as e 4iA . For 
reference, we record here the dispersive estimates that we will require in Section [U Note that 
these are specialized to the three-dimensional case. 

Proposition 2.9 (Weighted dispersion estimates). Let V : M 3 — > M be a potential satisfying 
HI, H2, and H3. Then 

lle^/llLK^^^IK^O*/!!^^), (2.23) 



/lk(R3) < -^-TTll(^)^7IUl( K 3), forpe [2,6]. (2.24) 



and more generally 

t 

Proof. The second statement (|2.24p follows by interpolating (|2.23p with the conserved L 2 norm. 
To prove (|2.23p . we note that 

W„itH f\\ llr>„— itAf 



7IU« = llfie-^nvil 



< 



e- itA STf\\ 



L6 

<- t \\(x)n*f\\ L 2. 



The first inequality is a result of the boundedness of discussed in Section 12.21 while the second 
is a standard dispersive estimate for the Schrodinger semigroup. Note that ||(x)fi*/||£2 ~ 

ll/ll^ + U^/HlUa. □ 

3. Multilinear harmonic analysis with potential 

In this section, we develop the multilinear analysis needed to build a space-time resonance 
theory in the inhomogeneous setting. We start in Section 13.11 with the basics, and proceed 
to prove in section 13.21 the analogues of Coifman-Meyer theorems for the distorted Fourier 
transform on M. d (d > 2) as stated in Theorem ll.il Then, in Section [3.31 we prove estimates for 
the commutators of the wave operator and position operators, such as [|x|, Here, we restrict 
ourselves to ffi 3 as we rely on Yajima's explicit description of the wave operator Q in [48] . We 
expect similar results to hold in other dimensions using the work in [491 [27] . These estimates 
will be instrumental in Section [3.41 where we attempt to understand the behavior of derivatives 
of the distribution M. that was discussed in Section 11.31 
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3.1. Definitions and first results. We start by considering pseudo-product operators of the 
form: 

T(f,g)(x):= [ [ m (6,6)/ tt (6)2%)e(*,6M*,6M6^6 (3.1) 
JR 3 JR 3 

Note that the case m = 1 corresponds (up to a constant factor) to the product of / and g. We 
say that the multiplier niQ satisfies Coifman-Meyer type bounds if the following homogeneous 
bounds hold for sufficiently many multi-indices a and /3: 

1^^^0(6,6)1 < cdai + i6ir (H+l/3|) (3.2) 

We would like to prove estimates of the form: 

ll r (/)3)llir'( K d) <m ,V \\f\\LP{R d )\\9\\Li{R d ) ( 3 -3) 

whenever f,g G <S(R d ) and p- = - + |. By duality, this is equivalent to proving that, for all 
/, ff , /» G S(R d ) and p, q, r satisfying 1 = 1 + 1 + 1 ( and = 1 - 1), 

A(f,g,h):= / T(f,g)(x)h{x)dx < \\f\\ifrBd)\\g\\ L i(R*)\\HLr(R*y (3-4) 
The left-hand side above can be written 

f T(f,g)hdx=[ [ m(6,6)/ S (6)<7 tt (6) / ^)e(x, 6)e(x, 6) #2 

JR d JR d JR d JR 3 

=■11 m(6,6)/H6)s%X^(6,6),^M6d6 

JR d JR d 

where we denoted by .M(6,6) £ <S'(M d ) the distribution determined by 

(-M(6,6),0) = f e(x,6)e(x,6)^ r V(^)^, </>GS(]R d ). (3.5) 

JR d 

One way to see that this is well-defined is to recall from ()2.9[) that 1 = OJ 7 ^ 1 , and hence 
by the boundedness of f2 on L 1 , J 7 " <^> G L 1 (R d ) whenever G 5(M d ). We remark also that if 
V = 0, then .M(6,6) = ^(' + 6 +6)- For convenience, we will abuse notation throughout and 
denote 

(.M(6,6),0> = / <K6)M(6,6,6)d6- 

As a result, we will write 
A(f,g,h)= [ T(f,g)hdx= [ [ m(6, 6)/*(6)<7 tt (6)^(6).M(6, 6, 6) #i dfc #3 (3.6) 

with the understanding explained above. 

We will also be interested in generalizations of (|3.6[) given by 

A(f,g,h)-.= [ [ m(6,6,6)/ fl (6)9 tt (6)^(6)^(6,6,6)d6^6d6- (3.7) 

■/R d JR 1 * 

where m is also a Coifman-Meyer multiplier in the three variables 6,6, 6 m the sense that 

1^6 "9^(6,6,6)1 < + 161 + |6ir (H+l/3|+l7l) . (3.8) 
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Our aim will be to prove the following two estimates on A(f,g,h): For p,q,r,p,q,r G (0, oo) 
satisfying | + | + £ = 1 and i + i + i = l + e, we have 

\Mf,9,h)\ <C aM ,v,e \\f\\L*{^)\\g\\LP{^)\\ h \\v-{^) (3-9) 
assuming the boundedness of 9\ on L q ,L p , and _/; and 

9,h)\ <C aM ,V,e H/llL<7(Md)|bllLP(R d )ll /l llL''(IR 3 ) + II / llL«(K d ) llfl 1 ll_.P(IR d ) II Hl* (M d ) (3.10) 

under no boundedness assumption on £R. This directly gives Theorem II .11 by duality. 

Remark 3.1. As will be evident from the proof, the second term on the right-hand side of (|3.10p 
can be further refined by projecting at least two of the functions f,g, h onto frequencies \D*\ < 1. 



Remark 3.2. Suppose that m is a Coifman-Meyer symbol as in (|3.8p . For each s > 0, we may 
define a symbol m s by 

m s (6.,6,£3) := "i(\/s6 ) \/s6,v / s6) • 

TTien m s is i£se// Coifman-Meyer with the same constants C a p~. We may therefore replace m 
by m s in either of the previous two results, and the resulting estimates will hold uniformly in 
s > 0. 

The remainder of this subsection is devoted to the proof of estimates (|3.9jl and (|3.10p . Before 
we begin, though, we will need the following maximal and square function estimates. 

Lemma 3.3. Suppose that W is an operator that is point-wise bounded by an LP -bounded positive 
operator, i.e. satisfying the point-wise bound 

\Wf(x)\ < CW\f\(x), for all f G L p (R d ), x G R d , 

for some positive operator W that is bounded on L p (R. d ) for 1 < p < oo. 

(a) Let tp G Co°(M. d ). For each n G M d , the operators 

/M- sup \We 2 ^^)f\ and f^ sup \We 2m ^ H^r)f\ 

iVi>iV 2 

are bounded on LP for all 1 < p < oo with a bound < (n) d . 

(b) For each n G M, d , the operators 

sup |e 2m — ^( — )/| and sup | e 2m ^ ^ ( }/ | ( 311 ) 

7V 6 2 Z iV Afi,Af 2 S2 z iV 2 

JVi>AT 2 



are bounded on L P (K ) for all 1 < p < oo wii/i a bound < (n) . 
(c) Lei U be any bounded operator on LP for some 1 < p < oo and suppose that {f n } C 
is a sequence of functions. Then 



IKEl^/nH^H^^IKEl/nl 2 ) 172 !!^ 
whenever the right-hand side is finite 
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(d) Moreover, if 4> is smooth and supported on an annulus, the operator 

sup \e 2 ^<pAf\ 2 ) l/2 (3-12) 

is bounded on LPi^ 1 ) for all 1 < p < oo with bound < (n) d '. 
Proof, (a) Since \Wf\ < W| /| , this follows directly from the fact that the operators 
/_^ SU p| e ^^)/| an d /-> sup \ e ^ iZ W^)f\ 

N ^ N 1 >N 2 iV 2 

are bounded on L p (M. d ) for 1 < p < oo. The latter follows from the point-wise inequality 

m^)f(x-y)\<(N\y\) d Mf(x), (3.13) 

where M f is the Hardy-Littlewood maximal function. The proof of this inequality is elementary, 
but we include it for completeness. Assuming without loss of generality that / > and picking 
a radial majorant for jj 7-1 ^) < X^<o ^Xb r (o) with J2j<o^ R j < °°> we can b° un d: 

1>(Z)f(* - V) ■= N d f F- 1 i,{Nz)f{x-y- z)dz = N d [ J c ~ 1, (/j(N(z — y))f{x — z)dz 

iV ./iBd .had 



<J>'W f{x-z)dz 

j<0 JB R 3 /N{y) 

V{N\y\) d i fix - z)dz < (N\y\) d Mf(x) if \y\ > 99Rj/N 



< 



(N\y\) d Mf(x). 



(b) Recall that by the intertwining property (|2.19p . m(D^)f = Qm(V)g where g = £l*f- Since 
Q* is bounded on LP for all 1 < p < oo, it is enough to show that the operator 

/ ' y SUp \n e ^^^)g\ 

is bounded on L p (M. d ) for all 1 < p < oo. But this follows from part (a) and Remark 12.61 

(c) Let U be a bounded operator on L p , for 1 < p < oo. By density, we may assume that {f n } 
is a finite sequence with n G {—L, ■ ■ ■ ,L}. Let {e n } be a sequence of i.i.d. random variables 
with P(e n = ±1) = 1/2 for n € {—L, . . . , L}. By Khinchine's inequality we have the point-wise 
equivalence 

(£ \Uf n \ 2 ) p/2 ~ E[| £ e n C//„n =E[\UJ2 tnfn\ p ]. 

n n n 
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d 



Here E denotes the expected value. Now, integrating both sides over R , we get using Fubini's 
theorem that 

<e[ ix;^/»i p ^~ii(Ei/»i a ) 1/a n^- 



(d) We start by noticing that the operator 

E l-Kin- V V7 
( sup |e— i>{ —)f\ 2 ) 1 ' 2 (3.14) 

is bounded from L p (R d ) to L p (]R d ). This follows by writing for any N\ > N%: 
\e * ^)f\ = \ e Nl ^( W )fN 2 \<{n) d Mf N2 

1 /2 

by part (a). The boundedness of (^2pj 2 \MfN 2 \ 2 ) is the Fefferman-Stein maximal inequality 
(cf. [121 143] ). The statement in (d) is now a direct consequence of part (c) with U = il. where 
0, is the positive operator majoring £1 as in Remark 12.61 □ 

3.2. Proof of Theorem ll.il With the maximal and square function estimates established, we 
are now prepared to prove the first main theorem. 

Proof of Theorem We will only consider the situation when m satisfies (|3.8p . The other 
case is similar (in fact easier). 

Step 1: decomposition of A. We start by decomposing f,g,h into Littlewood-Paley pieces with 
respect to the distorted Fourier transform as in (|2.20p . 

/= ^2 fa* 9= Yl 9N2 > h= Yl Hn 3- 

Ari62 z N 2 &2% A r 3 S2 z 

As a result, we get that 

Hf>9,h)= ^2 A (fNi,9N2,h N3 )- 

Ni,N 2 ,N 3 

By symmetry, we may assume without any loss of generality that N% < N2 < N±. Abusing 
notation somewhat, we will continue to denote the corresponding sum by A. Next, we let 
4> £ Cq°(M.) be given such that (f>4> = <f>. Define m Nl by the relation 

m«„&,6)«|w|w|)-^.(|,|,|). 

Then m Nl G C}°([—K/2, K/2] sd ) for some constant K depending on <f> with a uniformly bounded 
C 10rf+1 ° norrnj (independent of Ni, N2, N3). Expanding m Nl in a Fourier series, we can write, 



^Of course, this much regularity in m is not necessary, but it is sufficient for our purposes. 
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if (6,6,6) G [-K/2,K/2] 3d , 

^(6,6,6) = E « Ari K,n 2 ,n 3 )e¥K-5 1 +n 2 .6 + n3.6) 



m^ 1 



ni ,n 2 , 

where a^ 1 satisfies the bound: 

|a(ni,n 2 ,n 3 )| < m , (1 + |m| + H + \n 3 \)- 10d . (3.15) 

Consequently, 

A(f,g,h)= E E aJVl (ni,n 2 ,n 3 )5^X4 ( 3 - 16 ) 

N 3 <N 2 <Nx m,n 2 ,n 3 eZ d 

where 



-ni,n 2 ,n 3 
'N U N2,N 3 




e&*«*fUtij) (e^^te) 

.^(6,6,6)^6^6 (3.17) 

/l 1 ,n 1 (6)ff7 V2 , n2 ,7V 1 (6)^V 3 ,n3,JV 1 (6)-^(6,6,6) ^6 #2 
= // A 1 ,n 1 (6)A 2 ,n 2 ,Ar 1 (6)r(^iV3 1 n3,iV 1 )(6 ) 6)^1^2- (3-18) 

Here we have denoted 

^7Vi,ni - — e </ ATi > yN 2 ,n 2 ,Ni ■— c yjV 2 ' n N 3 ,n 3 ,N 1 -~ e '^3' 

and 

r(/iAT 3 ,„ 3 ,Ar 1 )(6,6) := (- A/i (6,6),^l r 3 in3i Ar 1 ) = / e(x, fi)e(x, b)h Ns ,n s ( x ) dx - 

J M. d 

Writing 

I. _ \ 1 ^ni,n 2 ,".3 _ \ ^ni,n 2 ,n 3 

n <N 2 ,n 3 ,N 1 — n N 3 ,n 3 ,N 1 , ^N lt N 2 ~ Z^/ ^N 1: N 2 ,N 3 

N 3 <N 2 N 3 <N 2 

we may further simplify (|3.16p to obtain 

\Hf,9,h)\< E E l«K,n 2 ,n 3 )p™ 3 | 

N 2 <Ni ni ,n 2 ,n 3 eZ d 

< m ,0 E a + M + Ki + Kir 1M - 10 E i s ™ 3 i- ( 3 - 19 ) 

ni,n 2 ,n 3 £l, d JV 2 <iVi 

As a result, we are reduced to proving the following estimate: 



| CT ni,n 2 ,n 3 | < , yd, yd, yd J \\J llpll«/llgll"llr for Part (i), 

'~ ( l} ( 2) ( 3) + for part (ii). ^ 
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Step 2: Proof of (|3.9p . We start by computing using Green's formulae: 

l£i| 2 r(£i,6) = / He(x,£ t i)e(x,£ t2 )h<N 2 ,n 3 ,N 1 { x ) dx 
JR d 

= / e(x,€i)H[e(x,Z 2 )h<N2 < n3,N 1 (x)]dx 

JR d 

= |6| 2 r(a,6)+ I e(x,Z 1 )e(x,&)H[h< N2>n3jNl (x)]dx (3.21) 

JR d 

- / e(x,€i)e(x,&)V(x)h<N 2 ,n 3 ,Ni(x)dx 

JR d 

+ 2 1 e(x,£i)Ve(x,&)V/i<Ar 2 ,n3,jVi(z) <te 

Thus, 

=// ^2/A, 1 ,„ 1 (6)l6| 2 4,n 2 ,7V 1 (6)r(/ i < J V 2 ,n 3 , J V 1 )(6,6)^1^2 (3.22) 

+ // |^/lr 1 ,„ 1 (ei)4,n 2 ,iv 1 (6)r(^<Ar 2> „3 I ivj(ei,e2)darf6 (3.23) 

+ // |^/i 1> „ 1 (ei)5k,n 2 ,AT 1 (6)r(y%iV 2! n3 > ivJfe,6)^ld6 (3.24) 
+ 2 /// ^/A 1 ,n 1 (^i)9l 2 , n2 ^ 1 (6)e(x,ei)Ve(x,6)V/ l < 7 v 2 , n 3, 7 v 1 (ix^i(i6. (3.25) 
We start with the contribution of ()3.22[) : Denoting by 



^^9N 2 ,n 2 , Nl = P jf)^9i with m := \£\ 2 (b(0 (3.26) 



and 



== ^ $7*fNum = & e^^±{-±-)J*f, with 0(0 := I^W)- (3-27) 



We now estimate 



£ rai< £ (^Y / / /! Vlini (Ci)4,n 2 ,AT 1 (e2)r(/ i <iv 2! n 3 ,iv 1 )(ei,e2)dei^ 2 

2 „„^ pjj 



N 2 <N! 



L Nl ,nS xS >3N2,n2,N 1 (x)(e KN i i/j(—)h(x))dx 



f „ .n 3 -D$ r)tt 

E^n/ME SU P \0»*,r»,»i\ 2 ) 1/7 ( SU P e 2m ^T^( — )Mx))dx 
jR d ^ - iv i."i ^ Ni>N 2 N!>N 2 ^2 



(ni) d (n2) d (n 3 ) d \\f\\ Lq{Rd) \\g\\ LP{Rd) \\h\\ Lr{Rd) . 
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The last inequality above follows from Lemma 13.31 The bound on (|3.23|) is similar: Denoting 



E rai< E [ [ /!v 1 n 1 (^A,n 2 ,JV 1 (6)r((-A + y)/ i <JV 2 ,n3,iV 1 )(Cl,6)^ld6 



N 2 <N r N 2 <N r 

N 3 <N 2 <N 1 v 1 



/jVi,m (x)gN 2 ,n 2 ,N 1 (a;)^jV3,n 3l JVi (:c)Gte 



(El/AT iril | 2 ) 1/2 ( SU P l^n^l) [ E ^ 



1/2 



JV 3 ,n3,iVil 



r/.r 



<( 



Ni 



(ni) (n 2 ) (n 3 ) 1 1 / 1 1 z,<? (Rrf) 1 1 1 1 (K^ ) 1 1 ^ 1 1 jl- (m^ ) , 



by Lemma 13.31 

We move on to studying the contribution of (|3.24p . Let s±, s 2 , S3 > satisfy s\ + s 2 + S3 = 2, 
si < d/g, S2 < d/p, and S3 < d/r. Then, 



E ®M\< E N ^ (^)5l 2 ,n 2 ,^ 1 (e2)r(^< J v 2 ,„3,J V J(6,6)^1^2 



N 2 <N! 



N 2 <N! 



E 



TV-, 



( Jv r"iv„„.)( JV 2"* 2 9»,™,iv 1 )(Jv 2 -«'»<K 2 ,„3,«,)i'W^ 



s/(E^r 2 *'i/ Nl ,„,i 2 ) 1/2 (E^ 2 " 

^ JVi jv 2 Nl ^ A 



< 



Ni>N 2 

. n 3 -pt ntt 
( sup N~ S3 e KN i ip{—)h)V{x)dx 

N 2 ,N! N 2 



\9N 2 ,n 2 ,Nj 



|2U/2 



<m»2»3> 



1'«C1B3>II5| 



< 



(ni) d (n2) d (n3) d ||/|| i9(Rd) || 5 || iP(Rd) ||/i|| Lr(Rti) 



where we denoted - = - = i = ^ and used Lemma 13.31 and Sobolev 

q q a ' p p a ' r r a ' 1 

embedding. We also remark here that in order to control (sup^r 2 ^ A jr 2 ~ S3 e KJV i tp(jj-)h) in 



terms of 



, we repeated the proof of part (a) of Lemma [3731 and bounded: 



|V| S , , V 



^)/i = i E ( 

JV 3 <JV 2 



NS ^N 2 



1\ 2 iV 3 



where ft := \^\ s (^(0 - V(2£)) G C °°(K 3 ) (recall |0(V/JV)/| < M/). 

Finally, we bound the contribution of (|3.25p under the assumption that 9\ in (|2.2ip is bounded 
on U>{ 
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N 2 <N 



E mB\< E ^r 2 Lu^s. 

N 2 <N ± JK 



,™2, 



dx 



N 3 <N 2 <Ni 

E 

N 3 <N 2 <N 1 

E * 

AT 3 <JV 2 
J62 N 

A r 2S2 z 
J,LG2 N 



A 2 ,n 2 ,Ni 



^2^3 



/ ^/Vi r)1 ( X )^5W2,n 2 ,iVi^A r 3,n3,A r i^ 



J 2 iVo 



7iV 2 ,n 2 , 9^^Af 3 ,n 3 , JiV 2 



r/.r 



l JN2tni ( x )^9N 2 ,n 2 ,JN 2 ^h N2 / Lin ,^ JN2 dx 



^ E 72T / ^pI/^^JEI^^^I 2 ) 172 ^!^^,^! 2 ) 172 ^. 

J,L£2 N ,/Rd ^ ' ^2 ^2 

Using the boundedness of 9^ along with Lemma |3.3|, we conclude from above that 

1 E 4223)1;$ E wMW^l/iiwiwi^iwu*. 



iV 2 <iVi J,L62 r 

This finishes the proof of part (i). 



(3.28) 



Step 3: proof of part (ii). Here we make no assumption on the boundedness of the Riesz trans- 
form 9\ associated to V. Instead, we assume that the potential V satisfies assumption H3*, 
which implies the L p (M. d ) boundedness of the operator *B defined in (|2,22|) . 
We split the analysis into three cases, depending on the size of Ni and N2: 
Case 1: N± > 1 and N2 > 1. In this case, we compute using Green's formula as in (|3.21|) . 
The only departure from the proof of part (i) is in bounding (13.25[) . for which we write: 



E <ra>|< E N l 2 [ l Nl , m ( X )^9N 2 ,n 2 , Nl V(h* N3 ^ Nl )dx 
l<N 3 <N 2 <Ni Jr3 



N 2 <Ni 
N 2 >1 



< 



E 



l<N 3 <N 2 <Ni 



N2N3 



where we denoted h* Ns ri3 Ni = /ijv 3 ,n 3l JVi wnen ^3 > 1 and h\ n3 Nl = h<i,n 3 ,iVi- The result now 
follows in this case exactly as in (I3.28|) . Case 2: N± > 1 and N2 < 1. In this case, we bound 
(with g = (Dt)P< ig ,h= (Dt)P< x h) 
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E rai< E ivr 2 / i Nltni (x)v 9N2 , n2 , Nl v(h N3 



N 2 <Ni 



l N 3 ,n 3 ,Ni 



N S <N 2 <1<N! 



N 3 <N 2 <1<N 1 



/at, r)1 ( 3; )^ r 2 5S 5A r 2,n2,AriA r 3 eS £h N3tn:j , Nl dx 



l 2 ) 1/2 <fe 



^ E ^ / IZ^ 1 , n J(E iV 2 e | , B^ 2 ,n 2 ,JV 1 | 2 ) 1/2 (E Ar 3" 2e |«^3,n3,iV 1 | 
JV X >1 ^ ' iV 2 N 3 

<l|su P |/ ||| L ,- sup ||( ^I^^^^^IIlpIKE^I^^^jvsI 2 ) 172 !!^ 

N - 7V ' ni K>1 ^ ^ 

where we have used Lemma 13.31 (part c)) and chose f* such that: 

lie 11 1 

— = and — I 1 = 1. 

r* f d q p f* 

Using Lemma |3, 31 and Sobolev embedding we arrive at the needed estimate: 
I E (E25])|<(ni) d (n 2 ) d (n3) d ||/|| L9 -|bllLH^IlL- 

N 2 <Nt 

Case 3: N% < 1. Here we refrain from doing the integration by parts after (|3.19[) and 
estimate: 



|^ni,n2,n 3 
I 1- 'JVijJVa 



En 2 r 
(aT) / ( N l\fNi,m\)( N 2 £ SU P l5JV2,n 2 ,JVi|)( SUp \h<N 2 ,n 3 ,Ni |) dx 
N 2 <Ni<l 1 3 Nl ^ N * N ^ N * 

%f (E ^llM.njWE^ 26 sup |5iV 2l n 2! iV 1 | 2 ) 1/2 ( sup I^^^Ddx 
7k3 ^ ^ JVi>JV a Wi>^V 2 



<(ni) d (n 2 ) d (n 3 ) d \\(Dt) e f< 6 4 L 4(Dt) £ g\\ Al -^ 
<(ni) d (n 2 ) d (n 3 } d \\f< 64 \\ L 4g\\ LP 



□ 



3.3. Commutators of position and wave operators. The purpose of this section is to prove 
the following theorem. 

Theorem 3.4. Assume that V satisfies assumption H3** as in Yajima [IE] (cf. Remark 2.7.) 
Then for any radial function a{x) satisfying |Va| < 1, the wave operator £1 satisfies: 



1 1 



.1 1. 



[a(x),Q] : L P (R 3 ) ->■ L q (R 3 ) for any 1 < p < q < oo m'tfi - - - < min(-, -); (3.29) 

p q 6 p 

In particular, [a(x),Q] is bounded from L p to L p+t whenever 1 < p < oo and e is sufficiently 
small. 

Remark 3.5. Typical examples for a(x) are \x\ and (x). In the latter case, the above theorem 
gives the endpoint case of Corollary 1.6 of [6], where Q. is proved to be bounded on (x)~P L P (M 3 ) 
forO<p<l. 
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Corollary 3.6. Assume that V satisfies assumption H3**. The wave operator Q is bounded 
on weighted spaces {x)~Plp(R 3 ) forO<P<l. 

Remark 3.7. Since [a(x),Cl\* = [f2*,a(x)], the above theorem is equivalent to the boundedness 
of [a(x), fi*] with the same range of exponents and same conditions on a. In fact, it is the latter 
statement that we prove below. 

Proof. By the remark above, the boundedness of [a(sc),0] in the (p,q) range specified in (13.291) 
would follow by interpolation once we show that 

[a(x), fl*} : L 1 (R 3 ) -> L p (R 3 ) (3.30) 

[a(x),Q*] ■ L q (R 3 ) -»• L°°(]R 3 ) (3.31) 

are bounded for any p £ (1,6/5) and any q > 6. We start with the first estimate (|3.30p . For 
this we need the expansion of the wave operator as 

fl*f = f- W x f + W 2 f - W 3 f + Lf (3.32) 

where W\, W2, W3 and L are bounded operators on IP for all 1 < p < 00 and will be described 
in more detail below. 

Step 1: Bounding [a(x), Wj]. Recall from [48} Equation (1.28)] that Wj takes the following form: 
Wjf= K j (ti,...,t j _ 1 ,r,u:)f(x + p)dt 1 ...dt j _ 1 drduj (3.33) 

J{0,ooy- 1 xIxY,J 

where S = S 2 (unit sphere in M 3 ); u = (uj\, . . . ,Uj) E £ J ; / = (—00, — 2ujj(x + t\0Ji + . . . + 
tj-iUj_i) is the range of integration of the variable r; y = y — 2(u)j ■ y)ujj is an isometry; and 
p = t\uj\ + . . . + tj-\Uj-\ — TUj. Furthermore, the kernels Kj (j = 1, 2, 3) are given by: 

K j (t 1 ,...,t j ,u>)= / e~ i ^-i^ tlSl/2 (si...s j )K j (s 1 uj 1 ,...,s j u} j )dsi...ds j (3.34) 

i[0,oo) n 

where 

j 

Kj(h, ...,kj) = i j (2iry 3j / 2 2- j Yl V(h - and k := 0. 

1=1 

For instance, 

Kliti,^) = c [ e- i ' lSl/2 siy(s 1 wi)(is 1 

J[0,oo) 

which can be written by integrating by parts (using the identity e lts l 2 = ^d s e tts / 2 ) as: 

Since, by Plancherel's theorem, K\ is also obviously bounded in L? Ul , we can write that: 

(tl) 2 |Xi(t!,a;i)| ^vPituuij+j/itum). (3.35) 
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where p € and p' £ Lf w . A similar calculation for Kj shows that if we define 

kj(h,. . . ,tj,ui,. . . := ^IJ(*i) 2 ^ \Kj(t 1 ,...,tj,U3 1 ,...,u}j)\, (3.36) 

then kj(t\, . . . , tj, u) can be bounded by a finite linear combination of products of the form: 

j 

l[pi(ti,u>i) (3.37) 

l=i 

where each pi(ti,ui) is either in L^ Ui or in L 2 u ■ Consequently, from now on we will assume 
without any loss of generality that: 

f \Kj (t u . . . , tj, u) | = nil (ti) ~ 2 kj (h, . . . , tj, u) q 
[kjih,. . . ,tj,u) <p l (ti,u l )...pj(tj,ujj); Pi(U,uji) e L^ u . +Ll i>0J .. 

Remark 3.8. It follows from the integration by parts argument leading to fj3.38f) . that the highest 
order term in K\ is exactly CV(0)/t 2 . This explains the failure of the boundedness of (I3.29|) 
when p = q. 

Changing variables in (|3.33|) by r i— > tj = — r — 2u)j(x + t\UJi + . . . + tj-iLJj-i), we get that: 
w jf( x ) = / K j (t 1 ,...,tj-i,-tj-a,u})f(x + p)dtdu} 

i[0,oo)J'xEJ 

where a = 2ujj(x+tiuji + . . .+tj-\Uj-i), t = (t±, . . . ,tj), and p = t\CJi + . . .+tjUJj. Consequently, 
we have the following expression for [a{x), Wj]: 

[a(x),Wj]= / [a{x)-a(x + p)]K j (ti,...,tj- 1 ,-tj-o-,u)f(x + p)dtdu (3.39) 

i[0,oo)J'xSi 

Now, notice that since a is radial and z i— )■ z — 2(z ■ uJj)uij is an isometry (reflection with 
respect to a plane), we have that a(x) = a [x — 2(x ■ ujj)ujj) which gives along with the bound 
on Va that: 

\a(x) — a(x + p)\ < \p + 2(x ■ = \p + 2{x ■ ujj)ujj — 2[(p + 2{x ■ ujj)ujj) ■ u)j]uij\ 

= \hu)i + . . . + tj-iUj-i - (tj + a)ujj\ < \t\\ + . . . + \tj-i\ + \tj + a\. 

We remark that it is important here to get a bound in terms of \tj + o~\ rather than \tj\, since 
the kernel Kj(t±, . . . , tj-\, — tj — cr, u;) in f|3.39|) does not decay in tj, but in tj + a. This leads to 
the following remark: 

Remark 3.9. From a technical point of view, if a was not assumed to be radial (e.g. a(x) = x), 
then \a(x) — a(x + p)\ can only be bounded by \p\ which would not give a bounded contribution 
from LP to L q even for the [a(x),W\]. In fact, an asymptotic expansion of K\ in (|3.34j) shows 
that its first order term is t~ 2 , for which the commutator [x,Q] does not seem to be bounded 
from IP — > L q whenever q > p, unless f is assumed to lie in some weighted space. 

This failure makes sense physically: 0, maps the initial data of the Schrodinger equation with 
operator H to its (free) scattering data at time +oo. From the quantum- classical correspondence, 
one can think of the solution at time zero as a particle coming from far away towards the 
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potential; it can be deflected for positive time and emerge from the interaction behaving like 
a free particle that started from a different direction. Thus Q has no obligation to preserve 
directions, and hence can only commute well with radial localization operators. 

As a result, we have: 

\[a(x), Wj]f\ < [ (|ti| + . . . + \tj-tl + \tj + a\) \Kjih, . . . , tj + a)\\f(x + p)\ dt du 

< f \ + 2 ' ' tJ* J ~> 2 J~ x Xf- kj(h, . . . ,tj + <r,ta)\f{x + p)\ dt du. 

J [0,00)3 X& {tl) Z ■ ■ ■ (tj-l) £ (tj + o) L 

(3.40) 

We will only prove the bounds (|3.30p and ()3.31|) for the operator 

<8if:= [ t— 2 + , . 2 k j (t 1 ,...,t j + a,u,)\f(x + p)\dtdw (3.41) 

as the others are similar. We assume without loss of generality that / > 0, and argue by duality. 
For any non-negative test function h, 

(<Sif,h) = [ / -—2 77 1 v2/ , , — rk j (ti,...,t j + a,u)f(x + p)h(x)dtdudx 

Jr* J[o,oo)jxSj \ti) ■ ■ ■ \tj-ir\tj +°~) 

< [f M £ w /•:•'• • /') (jiMtu^riur^ {tj+ar^t^p^+a^.ydtdudx} 1 ^ 



[O.oop'-ixRxSJ 



(3.42) 

j'-i 



x[/ h(xf f f(x + p)(]I(t l )^)\t j \ 2 dtd t vdx]W. (3.43) 

JRl J [0,ooy "^-xRxT.3 l=1 

We start by estimating (|3.42p . Recall that a = 2ujj(x + t\UJ\ + . . . + tj-iUj-i) and p = t\ui\ + 

. . . + tjUJj, so 



®M P =[ (l[pi(tl,"l) P (tl)- p ) f [ f(x + p) )^ [ v^+a^Ydxdtdu 

= [ (UPi(ti^i) p {*i)- p ) f [ f(v h } tj ^ 2P r- Pj (-t j + 2cj j .y,u; j ydydt ( ko 

J [0,00)3-1x1:3 \, J; JMJR3 {-t j + 2u j .y)P 



< 



U .\2-2p 

'jjRtJv* {-t j + 2u) j .y)P 
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since p < 2 and pi S + 1% u . As a result, we get that: 

r r \t-\ 2 - 2 P 

< ll/IUi sup / / 7 — UL - Pi (-^ + fc^.y, c^f dt, 



vm 3 J$?.. 7_! Jit.-IM (-U + 2ujj.y)P 



< 



duo 



if pj E Lt? >u) . + Lf. )UJ . and 1 < p < | thanks to Holder's inequality. Consequently, we have that 



» < n/iii( p . 



To estimate (I3.43p . we observe that 



sup/ /(x + /5) (IT^)" P ' J I*/ dtc ^ 

ceM 3 ^[0,oo)J- 1 xRxSi y ;=1 J 

/ / /(^ + h^i + . . . + tj-iujj-i + y) TT(fy)~ p dt\ . . . dtj-i du>i . . . dwj-i dy dx 



sup 



< 11 T, \L^ 



from which we can conclude 

< li h II .11 fiiW 



This shows that [a(x), Wj] : L 1 — >• L p for all p £ (1,6/5) and hence gives a satisfactory contri- 
bution to (l3T30j) . 

Next we show that this commutator is also bounded from L q — > L°° which verifies (13.31 H and 
concludes Step 1. As before, we only write the estimates for the contribution of (|3.4ip to (|3.40l) 
(the other terms are treated similarly). For this, we notice that for every x 6 R 3 , if 1 < q' < 6/5 
(or equivalently q > 6), then 

y d X 

®if(x) < [ [ f q (x + hcj! + ...+ t^) tjdtj duj J±- . . . 7 ii- dun... dLJj-l] 1 /* 

J[0,oo)ixEi \*l/ y \ t j-l) q 

7[0,oo),xS, (tl) q ' "' (tj-lY {tj+<j) ql 



£ WJWLl, 

by arguing exactly as when we estimated the contribution of (|3.42p above. 

Step 2: Bounding [a(x),L]. In order to prove the desired estimates for the integral operator 
L, we will need more information about its kernel L(x,y). These are included in the following 
lemma: 

Lemma 3.10 (|48|). The kernel L(x,y) of the integral operator L can be written as: 

LAx. v) — L_(x. v) 
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where L±(x,y) = Ya=o Z±,i{x, y) and Z± satisfy the following bounds: 

\Z±,i(x,y)\ < (M-lyl)- 3 ^)- 1 ^)- 1 fori = 0,2, 3. (3.44) 
\Z +tl (x,y) - Z_ tl (x,y)\ < (\x\ - \y\}~ 2 min ((x)^ 2 , (y)" 2 ) (3.45) 

Proof. All these estimate can be found in |48] (particularly in displays (3.10), (3.27), (3.28) and 
Propositions 3.2, 3.3, and 3.4). □ 

The operator [a(x),L] is an integral operator with kernel 

K(x,y) := [a(x) - a(y)]L(x,y). 

Since a is radial with bounded gradient \a(x) — a(y)\ = \a(\x\-^) — a(y)\ < ||x| — |y|| and hence 
by Lemma 13.101 we can bound: 

\K(x,y)\ <K 1 (x,y) + K 2 (x,y) 
where K\ and K 2 are the non-negative functions 

K 1 (x,y)=(\x\-\y\r 2 {x)- 1 {y)- 1 
K 2 {x,y) ={\x\ - |y|) _1 min((x)^ 2 , (y)~ 2 ) . 

Our aim is to show that the integral operators Tj : f(x) \-t J R3 Kj(x,y)f(y)dy are bounded 
from L 1 (M 3 ) -»• L P (R 3 ) and from L q (M. 3 ) to L°°(M 3 ) for 1 < p < 6/5 and 6 < q < oo. Since 
K(x, y) is symmetric, Tj is formally self-adjoint and it is enough to show that it is bounded from 
L 1 -> L p for all 1< p < 6/5. 

By Minkowski's inequality, this follows once we show that for j = 1,2: 

6 

sup||i^(-,y)|| LP (iR3) < 1 whenever 1 <p < -, (3.47) 
y a 

which is easily verified thanks to (|3.46|) . □ 

3.4. Derivatives of multilinear estimates. In this section, we develop an understanding of 
how to manage estimates involving the derivatives of A4 . As we discussed in Section 11.31 this 
understanding informs at a fundamental level the definition of space resonance in the distorted 
regime. Moreover, derivatives of M arise when one attempts to prove boundedness of the 
multilinear operators T and A (see (|3.ip and (|3.4p ) in weighted Sobolev spaces. In particular, 
we will see that in PDE applications, it is often important to be able to bound quantitates like 
xT(f,g) in L p , say. Since 

xT(f,g) = n*(p-\-iV^P)nT(f,g), 
and is bounded on L p , it suffices to study 

P^d^ [[ m(ei,6^3)/ ft (6)5 S (6)A4(6,6,e3)^i^2- (3.48) 



Naturally, this necessitates having a distributional interpretation for d^ 3 M. 



30 



P. GERMAIN, Z. HANI, AND S. WALSH 



To develop an intuition for what's going on, let us look at the case V = 0. There M = 
<5(£i + £2 + £3) and the above expression is equal to 

jr-i j_ j m(-6-6,6,6)^/'(-6-&)«7 l, (6)de2 

+ j d ?3 m(-£ 2 -£ 3 ,6,£3)/ 8 (-£3- 6)^(6)^2}- 

These are standard convolution integrals that one can estimate in LP in terms of Lebesgue norms 
of /, g and xf, for appropriate multipliers m. A moment's thought shows that this is possible 
due to the identity 

%<K6 + 6 + 6) = + 6 + &) 

which enables us to integrate by parts in £j. Consequently, we might hope that for V ^ 0, 

%-M(6,6,e 3 ) = d^M (6, 6, 6) + remainder, (3.49) 

where this remainder is a distribution of order 0. This is equivalent to saying that the vector 
field — preserves the order of the distribution M. . To test the validity chances of f|3.49[> , 
we inspect the case m = 1, for which ()3.49p would predict that 

£lxS}*(fg) = (Qx(Tf)g + remainder. 

The relation 

n x n*(f 9 ) = xfg+[n,x]n*(f g ) = (n x n* f) g + ([n,x]n* f) g + [n,x]n*(f g ), 

shows that the validity of the previous line is tied to the boundedness of the commutator [x, O] 
on LP. Remark 13.91 tells us, however, that [x,Q] cannot be bounded since x is not radial. It is 
then more reasonable to expect that 

Cl\x\Cl*(fg) = (fi I x\ Q* f)g + remainder, 

with the remainder bounded in LP in terms of Lebesgue norms of / and d^- In other words, we 
seek to relate |(9^ 3 |A^ := |V^ 3 |A^ and ^^.M = IV^J-M, rather than something of the form of 
(13.491) . This strategy turns out to be the right one, and we are able to exhibit precisely such a 
relation, both qualitatively in (|3.57p . and quantitatively in Theorem 13.131 

Remark 3.11. Ln what follows, we will use the summation convention that repeated indices are 
summed over from 1 to 3. For instance, \V X \ = R a d x a where R a = is the Euclidean Riesz 
transform. 

3.4.1. Understanding d^M . Let us start by looking at (|3,48p with m = 1. By duality, integrating 
(I3.48P against a test function h, we get 




#2 (3-50) 



It seems important for our future work (and much cleaner for our current application) that estimates on this 
remainder do not involve weights falling on / or g. 



RESONANCES IN THE PRESENCE OF A POTENTIAL 31 

Prom the identity x a f = ^(J 7 " id^a)^)Qf we have (dropping for convenience the constant 
-i) 

- (I330D = / f(x)g(x)n\x\n*(n^-n*)h(x)dx 

./lR3 \X\ 

\x\f(x)g(x)(n^-n*)h(x)dx+ [ f(x)g(x)[n,\x\]n*(n^n*)h(x)dx 

M 7r3 \x\ 

ol r ol 

{n\x\n*)f(x) g (x)(n^-n*)h(x)dx+ / ([\x\,n]n*f(x)) g (x)(n- ri n*)h(x)dx 



x Jm3 x 



f(x)g(x)\\x\,n]n* (n—n*)h(x) dx. 

\x\ 



Define the operators 



x 



n a :=n—n*, £ = [\x\,n]w 



and their Fourier space manifestations 

R a ■= PlZ a P~ X , E = P£p~ l . (3.51) 

so that J-^lZ a f = R a J^ f and similarly for £ and E. Notice that since J-^fi = J 7 , R a is nothing 
but the Euclidean Riesz transform i^- As a result, we may write 



(1330]) = / TZ f3 (nx/3n*)f(x)g(x)n a h(x)dx 

+ f £f(x)g(x)TZ a h(x)dx- [ f(x)g(x)£TZ a h(x) dx (3.52) 



R l a € ? & ^ (6 ) Rg h* fa)M (6 , 6 , £s Ri d£ 2 #3 (3.53) 
+ / ^/ S (6)5 tt (6)^^fe)-M(ei, 6,6Ri« 3 (3.54) 

iR3 

- I fKZi)9 i (b)ERght&)M(!i 1 ,b,&)dZ 1 dt; 2 d{; 3 , (3.55) 
In total, this gives the following formal expression for the distribution d^Ai: 

d^M =R%R^d^M - RgE^M + RgE^M or, (3.56) 

\V i3 \M =|V a |M - E^M + E* i3 M (3.57) 

where we use the notation Ee t to indicate that E is being applied in the £1 variable fixing 6,6- 
Recall that R a is self-adjoint and commutes with differentiation. 

Remark 3.12. All the entities in (|3.56[) and ()3.57|) can be interpreted as distributions on M 9 
as follows: Clearly, one can make sense of them when applied to test functions in 5(IR 9 ) that 
are linear combinations of tens ored functions (i.e. of the form </>i(£i)<^2 (6) </>3 {&))■ In fact this 
is explicitly done in ()3.53[) . (|3.54p . and (|3.55[) and the physical space manifestation in ()3.52[) . 
Since this gives a continuous linear functional on a dense subset ofS(M. 9 ), the relations ([3.56P 
and (|3.57p are well-defined in 5'(M 9 ) by density. 
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3.4.2. Multilinear estimates. After this qualitative formulation, we are now ready to obtain 
quantitative estimates on (|3.48p in U . By duality, integrating (|3.48|> against a test function h, 
we get: 



6, 6)/ tt (6)5 fl (6)%^(6)^(6, 6,6)^1^2^3, 



(3.58) 



which is equal to 



■/?? 



(6,6, 6)/ fl (6k%)^)%.M(6, 6, 6Ri <% 2 



+ / %^(ei,e 2 ,6)/"(ei)5*(6)^(6)^(ei,6,6)^i^2^3 



m 



(6, 6, 6)/ tt (6)s s (6)^(6)l -M(6, 6, 6R1 d& d& 



[™(6,6,6)/"(6)<?"(6)^(6)l -^(6,6,6)^6^6 



+ / 



m(6, 6, 6)/ tt (6)<?'(6)^(6)J -M(6, 6,6)^6 dfc 
+ / %m(6,6,6)/ s (6)5 tt (6)/i tt (6)^(6,6,6)^6rf6d6- 

With some rearranging, this becomes 



(1338]) 



%W«(6)^(6) -^ Q 3< (^™/ 8 (6)^(6))] <? S (6)M(6, 6, 6)^6^6 

(3.59) 



(6, 6, 6)^/ s (6)5 tt (6)^(6)] M(6,6,6Rid6^6 



RgE^ [m(6,6,6)/ s (6)5 tt (6)^(6)l -^(6,6,6)^6^6 



+ / E e ,R 



^(6,6,6)/ s (6)5 tt (6)^(6) ^(6,6,6)^1^2^3- 

This suggests looking at the following operators 

Ai :(/,<?,/*) ^ / rn(6,6,6)/ s (6)5 tt (6)^(6)^3<^(6,6,6) 

Ai:(f,g,h)^ [ m(6, 6,6)^(6)^(6)^(6)^^3^(6, 6,6) 

A 3 :(/,<?» ^ / m(6,6,6)/ s (6)5»(6)^(6)^3^ 3 ^(6,6,6) 



(3.60) 
(3.61) 
(3.62) 

(3.63a) 
(3.63b) 
(3.63c) 



Theorem 3.13. Suppose that V is a real potential satisfying HI, H2, H3*, and H3**. Let 

m be a symbol of the form 



m(6,6,6) = (161 + 161 + l6irV)(6,6,6), 



(3.64) 
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for a Coifman-Meyer symbol tjiq and some 5 > 0. For any q,p,r, q,p,r in (l,oo) such that 



111 111 

— I h - = 1 and - + - + -> 1, 

p q r p q r 



we have for all i = 1,2,3 



|Ai(/,5^)l<ll/llpll3l| g ||^l|r + ||/||p||5iy%- 

Proof. We only bound Ai, the argument for A2 and A3 being similar. We first split /, g and h 
into their distorted Littlewood-Paley pieces and write: 



A r i,A r 2,A r 3 62 z 



R l R l 



(6,6, 6) A, (6)<4 (6)4, (6)1 -M(6, 6, 6)^6 #2 d& 



"N 3 



We only present the calculation for the case when JVi > N 2 > N%. The other cases can be 
treated analogously. Expanding thq as before, we get 



™o(6,6,6M^¥(J^(J^) 



a Nl ( ni ,n 2 ,n 3 )e™i 



("■l-?l+".2-^2+n3-C3) 



"1,^2,7136 



where \a Nl (n±, n 2 , n 3 )\ < (1 + \n± \ + |tt-2 | + 1^3 1) 100 uniformly in N±. Consequently, the contri- 
bution of the case N\ > N 2 > N3 is bounded by 



5^ (l + |ni| + |n 2 | + 



n-3 1 



-100 



ni,n2,n 3 
X 



* — /" « 27rin 1 -g 1 r u 27ri n 2 'S 2 u 27rin 3 g3 

£ / <[e^^(6)^A 1 (6)](e^^^ 2 (6)) J Rf 3 (e^^4^(6)) 



It therefore suffices to prove that 



N2<N\ 



< 



^(6,6,6)^6^6^6 



(<ni><7z 2 >(n 3 >) 3 |mMl5l|L 9 



(3.65) 



Here we have denoted, following the previous convention, 



9N 2 ,n 2 ,N 1 



-1 2iri a 



h< N2 ,n 3 M x ) ■= ^ e^^4 7V2 (6). 

We divide the above sum into three pieces: first Ni,N 2 > 1, second Ni,N 2 < 1 and finally 
N 2 < 1 < N±. In fact, the last piece can be treated by combining in an obvious fashion the 
arguments for the first two, so we shall skip it here. For the first piece, we use (an elementary 
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consequence of) Lemma f3T3l to obtain 



l<iV 2 <JVi 



TZ 13 & /jVi,ni(aj) 9N 2 ,n2,N 1 (x)'}l a h<N 2 ,n3,N 1 (x) dx 

i$ Yl ^ S \\fNi,ni\\p\\9tf a ,na,Nx\\q\\h<N a ,n3,Ni\\r < ("-2) ("3) ) 3 II I\\lp MIli \\ h\\ L r . 

N 2 <Ni 

Turning to the second piece, let Q be denned by -k = 1 — i — A < i. Then from Bernstein's 
inequality we have that 

9N 2 ,n 2 ,Ni( x ) 'R> a h<N 2m> N 1 (x) dx 

N 2 <N 1 <1 •' R9 

~ ll/A r i,nillpllfi , A r 2,n2,A r illQll /j <A r 2,n3,A r lllr 
N 2 <Ni<l 

3(1 — — ) 

~ X/ ll/jVl.nillp-^2 9 Q ll5 , Af2,n2,Afillgll^<A r 2,ri3,A r i II? 
N 2 <Ni<l 

<((ni)(n 2 )(n 3 )nf\\ L 4g\\ Lq \\h\\ Lr . 
This completes the proof. □ 

4. Small-data scattering for a quadratic NLS 

4.1. The a priori estimate. The proof of Theorem 11.21 essentially consists of an a priori 
estimate in the space X given by the norm 

\\f\\ x -.= \\ x n*f(t, x )\\ LrL 2 x + \\f\\ LrH i ~ \\\x\n*f(t,x)\\ L?Li + \\n*f\\ Lrm , (4.i) 

thanks to assumption H3*. Due to the intertwining property of and (distorted) Plancherel 
identity, the X norm is equivalent to the following: 



Observe that Duhamel's formula for equation (jl.3j) can be written 

f = u -iB(f,f), 

where 

W, /)]*(*, £) := |y /" e-^ K '' ,,C) F(fl^)F(«,C)^(e^,C)^dCds, (4-2) 

with 7/, C) = |£| 2 + \rj\ 2 + |C| 2 . The key point is then to prove the following result. 

Proposition 4.1 (A priori bound). If f solves f = Uq — iB(f, f) and belongs to X with a norm 
less than a fixed constant e, then 

\\B(f,f)\\ x <\\f\\ 2 x 

The above proposition entails the inequality \\f\\x i$ H^ollx + ||/||x> which gives an a priori 
control on / in X. This leads to the main theorem by a continuous induction argument. In 
particular, thanks to the L 2 -convergence of the integral in (|4.2p . which follows from (|2.23p . we 
also get that / converges in L 2 as t — > 00 to some (j)^. This is equivalent to the scattering of 
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u(t) to a linear solution e ttH (j)^. By the boundedness of O on L 2 , one also obtains scattering to 
a free solution e^cj)'^ for some cj)'^ € L 2 . 

4.2. Decomposition of B(f, /). We will first decompose B into high and low frequencies: call I 
the piece corresponding to | (£, r], Q | > 1/ \/s, and II the piece corresponding to | (£, 77, £) | < 1/ -y/s. 
We then perform a normal form transformation on the high frequency piece (which is localized 
away from the time resonant set). Finally, we estimate separately the resulting terms Ii, I2, I3, 
I 4 

To be more specific, let x be a smooth cutoff function satisfying 



xi^^g-i if |£| 2 + + | C |2 > 4 



For any s > 0, we further define \s = x{y/s')i observe that the support of Xs will be a subset of 
^2/ v / i(^)- Using x to partition the frequency domain, B(f,f) becomes 

£?(/,/)» =1 + 11 

with 

I := J If (1 ~ Xs)e- is *TKv)P(OM& V, dr, dQ ds (4.3) 

II := f ff Xse- lS(f, p(v)TH()M(Z, V, dr, dQ ds. (4.4) 



where 



We now perform a normal form transformation on I (integration by parts using the identity 
-T$ d s e ~ isrp = e"^)- This gives 

I = Ii + I 2 + I 3 + I 4 

li := II ^^e-^JKvWOMdrjdC (4.5) 
h ■= ~ fj f ^e-^p(ri)p(()Md V d(ds (4.6) 

1 3 := fj f l -^e-^dJKri)7KQMdi 1 dQds (4.7) 

14 := fj f l -^e-^7Kil)ds7KQMdi 1 dQds (4.8) 

Here we have suppressed some dependencies in the interest of readability; we will continue using 
this convention throughout this section. 

Finally, it will often be necessary to distinguish which of the frequencies is largest. This will 
be the role of the cutoff functions = rj, Q, j = 1,2,3 that we now define. These 

functions are chosen to be smooth, and homogeneous of degree outside of the ball of radius 
1/10. Furthermore, they add up to 1 away from (0,0,0): 

$(1) + $(2) + $(3) = 1; for aU ^ Q with |£|2 + ^2 + |^|2 > 1 (- 4 Qa ^ 
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and are such that 



OnBitornSuppSW |£|>M + |CI 

OnJ3i(0) c nSupp*W, H>|e| + |CI (4-9b) 
OnJ3i(0) c nSupp*( 3 ), |CI>I£I + M- 

In line with our convention, denote 3>i = for each s > 0. 

4.3. A priori estimates: proof of Proposition [4TTI In this section, we prove Proposition ^. 11 
We begin with two simple observations. First, for any g 6 X, 

\\g\\L 2 (R3) ^ H^*slU 2 (]R3) ^ \\xtt*g\\ L 2( B cj + ||^*5||l2( Bi ) 

< \\xQ*g\\ L 2 (M 3 ) + ||0*5|| L 6( K 3) < ||X0*5|| L 2 (]K 3) + ||VO*0|| L 2 

< ||xO*5r|| L 2 (M 3) + \\D^g\\ L 2 



Thus, to control the full X norm, we need only prove \\xn*B(f,f)\\ L 2, \\D$B(f, f)\\ L 2 < \\f\\ 2 x . 
Second, for 6/5 < p < 6, we can bound the L^°Lx norm of a function in terms of the X norm 
by noticing that 

f IK^*/IIl-l2 if 6/5 <p < 2 1 

L°°Ll ^ II^/IIl-l^ S II nil / lie *n/i|i-e ^fo^^^fi ? < (4.10) 



\\ D *f\\ 6 L?Ll\\f\\ l LrLl if2<P<6 



Here, the estimate for 6/5<p<2is simply from Holder's inequality, and for 2 < p < 6 we 
interpolate between l? x and L^, controlling the latter with the (Euclidean) Sobolev inequality. 

Even at this stage, it is clear that our proof will require considering an inordinate number 
of terms: each of (|4,5p - (|4,8p may be carved into multiple pieces using the cutoff functions. 
Indeed, we will compound the problem by, at certain points, differentiating various terms as well. 
Thankfully, the vast majority of these quantities can be estimated using slight variations on a 
few basic strategies. With that in mind, we give a full treatment of each class of representative 
term, as well as the most problematic ones. In hopes of keeping our exposition to a readable 
length, however, when a quantity can be controlled using a straightforward modification of an 
already stated argument, we will only provide a sketch. 



Proof of Proposition pT"/} By the preceding remarks, we know that we must control B(f,f) in 
the weighted L? and norms in terms of Let us estimate each of the pieces Ii, . . . , I4, II 

in these spaces one at a time. 



4.3.1. Weighted L estimate of Ii for t > 1. We begin by controlling the long time weighted 
norm Hxfrj 7 " (|4.5j) ||£2 ~ ||ck (|4.5j) ||£2. Differentiating Ii in £, and using (|3.56j) . yields (sticking 
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to the convention that repeated indices are summed) 

agii =U ^mle-^JKnWKQMdvdt (4.n a ) 

* m 2 e -it4> ft^ ft (QRpRfMdr] d( (4.11b) 

l^e^d^MvWOtf^MdridC (4.11c) 

+ / / ^^e-^TKv)TKO [R%E V + B^] Md V d(. (4.11d) 



where 



m < = Vr + (1 - Xt)-T7r - K 1 ~ Xt) 



t(j) t(j)^ \J~^> 

{dr,PX)t d pcf) d /3(f) 

+ xt)—ri= - I 1 - xt)—r 



t(j) t(f)y/4> sf§ 

satisfy uniform Coifman-Meyer bounds and are supported outside of a ball of approximate radius 
1/y/t. Notice that (l47TT5|> - <|47TTd|> arise when we use (|3.56p to evaluate d^Ai, and then perform 
an integration by parts in rj. To simplify the notation, we omit from now on the superscripts a 
and (3 when there is no chance for ambiguity. 

We start our analysis with (I4.11ap . For any 5 > sufficiently small, it can be bounded by 

= \\tff ^mle-^^mTKOMdrjdCh, 

^tOle^A-Vllxalle^/IUe + He^A-Vll^lle^/lUe] 

<t\\f\\LS/AW tH f\\L* 

<\\ftx- 

Here we have used Proposition 12.81 (|4. l(Jj) . and expressed 

in order to apply Theorem 11.11 with the Coifman-Meyer symbols rf 1 j\f^>. This last step relies on 
the IP boundedness of the operator /\D^\, which follows from (|2.19p and the boundedness 
of the Euclidean Riesz transform. 

The term (|4.11bj) can be bounded similarly, so we turn directly to (]4. llcH . We need to 

decompose this term further, with the help of the cutoff functions x, 1 — x, &[ 2 \ ■ 

This gives six elementary pieces, of which we will only present the details for the ones involving 
$| 2) ; those with and $J 3) can be managed similarly, with the small modification of a duality 
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argument for . Thus we consider 

l^X^e-^dJi^TKOR^Md^dC (4.12a) 



+ 



/ [ -X^e-^fKrifKORffiMd-ndC (4.12b) 



To estimate (|4.12a|) . we use Theorem 13.131 and Proposition 12.81 for any 6 sufficiently small, we 
have 

||(|I32aJ)b 2 = || If (1 ~g )|7?|2 X^ 2) e-^^dJ^JKOR^RtMdrjda^ 

< \\e- itH A; 2 nxn*f\\ L s\\e- itH f\\ L e + ||e-^A t - 2 fix07IU3||e-^/|| L6 -, 

<t^\\n x n*f\\ L ^\\f\\ x <-^\\/f x . 

To estimate (|4.12bj) , also rely on Theorem 13.131 to get 

(1-xOM/., _ v u(2) -#ii, liUYfit 



\^m\\v = \\ II y -^^{i-xW ] e-^^d v fKv)fK0R^RtMd^dC\\ L 2 



< He-^A^nxfiVll^lle-^/IU- + ||e-^Ar 1 ^/||L3||e-^/|| i 6- 

^ll/lk < £7/8 



<||^*/|| i2 _||/|| x <_||/|||. 



The remaining term (|4.11dj) is treated the same way. 

4.3.2. Weighted L 2 estimate of I\ for < t < 1. Here the general idea is to repeat the argu- 
ments of the previous subsection, but systematically replace the dispersion estimates with the 
(distorted) Sobolev inequality of Proposition 12.81 (b). For instance, 

\\^m\\ L 2 = \\tjj K m l e -^^p {r]) p {0 M dvdCh2 

< t [\\e ttH A- 1 f\\„\\e itH f\\ L * + \\e UH A" 1 f\\ L3 s \\e itH f\\ L e] 

<t\\f\\ L5/ 4f\\ H} <t\\f\\ x 



Similarly, to estimate (|4.11b|) . we employ the cutoff functions. To bound a representative term 
of the form (|4.12ap . we note that 

\\@mh> = II // (1 ~$ )M2 *^ 2 ^-^T^dJ'WiOR^MdrjdQ 

< ||e-^A,- 2 0x^/llL3||e-^/|| L6 + \\e- itH At 2 nxn*f\\ L 4e- itH f\\ L6 -s 



<t^n x n*f\\ L2 \\f\\ H} <t^\\f\\ 



All in all, this gives the desired bound ||9^Ii||2 < \\f\\x f° r * < 1- 
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4.3.3. Weighted L 2 estimate 0/I2 for t > 1. First, we need to split the time integral giving I2 
into two pieces: 

ft 



So Si 



I2 + I|- 



The term 1 3, will be dealt with when we consider 1 2 for t < 1; thus for now, we only consider 
the term l\. Differentiating it in £ gives 

d^ll = j J J ^smle-^7Kn)7KQMdr,dQ (4.13a) 

+ J ff ^sm 2 s e- is *TKv)TK0RfatMdr,d( (4.13b) 

+ f ff mle- is %eJt( V )Tt(QRPR%Md V d( (4.13c) 

+ J JI ^~ iS ^TKv)P(0 [R^E V + B^R*] MdridQ, (4.13d) 

where V refers to the gradient with respect to (£,77, £). Here 

1 fl a^chvxM ^ . (e^,o-(vx) s r „ , 

V S( P J V S( P 
™. := { ^ J + » ^ v^V*' 

m s := (^^C)-(Vx), 

2y / SZ0 

are supported on annuli with inner and outer radii approximately 1/y/s, and satisfy uniform 
Coifman-Meyer bounds. We begin with (|4.13ap which can be controlled with the help of Theo- 
rem [Tj] and Bernstein inequalities: for 6/5 < p < 4/3 and 5 > sufficiently small 

IKSMdU* < 5 ^V^IIe^P^/HLslle^/IUeda 

+ j* ^sV sH P <l/v -J\\ L 4e lsH f\\ L .- S ds 



? \\ x ^' 14 ^ 



< / V~ss*--?\\f\\ LP ^ds 



The term (|4.13bp admits a similar treatment, thus we skip it and focus on (|4.13cp . Due to 
the presence of the (Euclidean) Riesz transforms, we will need to use Theorem 13.131 This is 
permissible because, for t > 1, enjoys the improved Coifman-Meyer structure of (|3.64p . 
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Thus, for t > 1, we may select 5 > small enough so that 



ds 



\\&m\\ L *<J II // m^-^/tt^O^TW^dCll^ ^ 

< 5 f [||e^P <1/Vi OxO*/||L3||e isif /|k6 + ||e^P <1/v/i ^*/||L3||e lsH /|| i6 -, 

</ ^jlinxjr/ii^J^d. 
< II/I&- 

Finally, (|4.13dp can be estimated similarly, thus we skip it. 

4.3.4. Weighted L 2 estimate 0/I2 for < t < 1. Once again, it suffices to follow the argument of 
the previous subsection merely replacing the dispersive estimates with appeals to the (distorted) 
Sobolev embedding inequalities. 

4.3.5. Weighted L 2 estimate 0/I3 for t > 1. First let us point out that, due to symmetry, the 
argument we give for I3 applies to I4 as well. The main idea in both cases is that we can use 
the equation satisfied by u to replace the quadratic nonlinearity with a cubic one. Indeed, this 
is immediately apparent since, by the definition of the profile / and the equation satisfied by u, 

dJ = d s {e- isH u) = e- %sH u 2 . (4.15) 

Using (I4.15D will give us additional decay, but it presents a technical nuisance in the form of an 
asymmetry in the integrand. We start by splitting the integral defining I3 into 

The estimate of I3 follows from the argument we give for I3 with < t < 1 in the following 
subsection; so let us focus on I3. Using (|3.56p . we see that differentiating it in £ gives 



€ I§ = -J JJ ^mle-^dsp^fHOMdvdCds (4.16a) 
+ fff ^m 2 s e- i ^d s p(v)TK()R(R'iMdr 1 d(ds (4.16b) 



+ 111 ^^e-^d s JKri)dJ~KOR%RtMd V d(ds (4.16c) 



where in this case 



+ f J f i^e-"*«>/»(>?)/»(0[S? B < + S 5 flf]jK dr, d( ds (4.16d) 



m * := — TH + ( X ~ Ks)——h- + (1 - Xs)^ 
satisfy uniform Coifman- Meyer bounds and are supported outside of a ball of radius ~ l/y/s. 
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First consider (l4.16aH . Introducing the functions and arguing as before, we wish to 

establish the L| boundedness of terms of the form 



sm s 



$®e- ia +d 8 ft(T,)ft(0MdridCd8 ) j = 1,2,3. (4.17) 



(2) 

We simply deal with the term involving & s , since it is possible to treat the other ones in a 
similar way. 

I|(S3Z|IIl|</^II // ^Pe- is *dJ-Kv)P(0Mdr,dC\\ L 2 ds 

< 5 [* s\\Aj 1 e- i ' H dJ\\ I j i (\\e i ' H f\\ Ifi + \\e lsH f\\ L es) ds 



<^ s ( s V4|| e -^ s /|| L2 )(-^||/|| x )^ 

< j\( S V^\\f\&)(-^\\f\\ x ) ds < \\ff x . 

Here we have used the fact that 

\\e- isH d s f\\ L * = \\u 2 \\ L 2 = \\e isH ff Li < s- 3 / 2 \\f\\ 2 x . 
The term (|4.16bp is almost identical, thus we immediately consider ()4. 16c|) . It can be split as 

M =E/' II // " xH j) e-^d s TKr 1 )dJ~KC)R^Md V dC\\ q ds (4.18a) 



j'=i 



+ E /* I' // ^^X^e-^dsfKnWKO^M dV d(\\q ds, (4.18b) 

We single out one of these terms, (|4.18bp for j = 2, the other being similar. It can be estimated 
as follows: 

||(M) J=2 || 2 <f\\JJ ^ ~ *0 HICI *^ 2)e ~^ dsTK^AJKOR^RtMd^dChi ds 



< 



{Aj'e-^dsf^ + WK'e^ 11 d a f\\ 3 _ 5 



< | >A[||« 2 ||3 + ||« 2 || 3 _i] ||nxn*/ll2tfa 



4.3.6. Weighted I? estimate of I3 /or < t < 1. Again, in light of the distorted Sobolev 
embedding theorem, the small time bounds are essentially a corollary of the long time bounds. 
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4.3.7. ifjj 1 estimates of I. Let us consider one at a time the terms Ii to I4. Multiplying Ii by £, 
partitioning in frequency, and then taking the L| norms gives 

H^Iilljri < II // ^^e-^JK^TKO^M d V dC\\ L 2 (4.19) 

+ II JJ t^^e-^TKvWO^MdrtdCh* (4.20) 

+ II // t^^e- lt ^)TK(M 3) Md v d(\\ LV (4.21) 

The procedure for each of these is similar, so let us consider only the second one. For small time 
we may estimate by Theorem 11.11 and the Sobolev embedding theorem: if t < 1, 

< II // iH^^e^^K^TKO^Mdr^daq 

< WK^fh'Qle^fWv + \\e itH f\\ L 3-s) 

<Vt\\f\\ 2 Hl <\\f\\ 2 x . 

The long time control is found in the same way only using dispersion: if t > 1, 

< W^^fUed^fhs + \\e itH f\\ L3 - S ) 
<Vt\\e itH f\\ L6 (\\e itH f\\ L s + \\e« B f\\ I t-,) 



< Vt-\\f\\ 2 x . 



1 
7 

Next, multiplying 1 2 by £, and then taking the L| norm, we arrive at 

ll^ rl i 2 ||^ ~ ^ I' // ^ C U e-^7»(^)7»(C)^w ^ czcil^i ^ 

r m.e-*+p(ri)7l(QMd7idC\\ L * da, 







where m s is a Coifman-Meyer symbol supported on an annulus of size ~ Vv^- For the small 
time control we use the (distorted) Sobolev embedding: if t < 1 



£ f ^ r V sH f\\ L A\e tsH f\\ L sds + r^lle^/II^He^/IUods 
» Jo V s Jo V s 



1 



< I ^||/|| 2 m(is < N 



'3 



The long time bounds are the proved in the same way with the dispersive estimates in place of 
the Sobolev inequality. 

Finally, terms I3 and I4 are equivalent by symmetry, so let us only discuss the former. We 
being by splitting the time integral as 
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For Ig, we partition the frequency space using the cutoff functions <&^\ j = 1, 2, 3. The resulting 
pieces are treated similarly, so we consider only the representative term with j = 2. In that 
case, we argue as follows: if t > 1, 

ll^-^llflj JJ^l-Xs^e-^dsf^fHOMdvdaqds 

<J'V~S [\\e isH d s f\\ L s\\e isH f\\ L e + \\e isH d s f\\ L s \\e iaH f\\ L e-s] ds 
< J [\/s||« 2 ||i3||it|| L 6 + Vs||n 2 || L3 ||u| 



L 6-<sJ ds 

1 1 1 



< / ^||/||3 _d S <||/|| 3 x 



X 



For, I3, we repeat the argument above substituting as usual the Sobolev embedding theorem for 
the dispersive estimates. Thus 



< 



/ [aAII^IIlsII^IIl 6 + v / s||« 2 IIl 3 II w IIl6-<5] ds 

Jo 

V^ll/llii ds<\\f\\ x . 

tt 

Notice that the same argument works to bound \\J$ 1 13 1 1 ^1 for t < 1. 

4.3.8. Weighted 1? estimate of II for t > 1. In contrast to the normal forms method used to 
control I, when estimating II we only want to use the fact that the integrand is supported 
compactly due to the presence of the Xs factor. Recall that 



n = J JJ Xse- ls *P(v)fKC)Md V d(ds, 

We split it as 

II = / . . + / . . =: II 1 +II 2 . 



We start with II 2 ; the term II 1 will be estimated along with the case t < 1 in the following 
subsection. Applying and using (|3.56p to evaluate d^M. gives 



-t 

d^II 2 = I 1 1 yfs~mle- is *p{ri)p(C)MdridCds (4.22a) 

^- sm 2 e -is^> p ^ p (Q R P R *J4 drj 4£ ds (4.22b) 

Xse-^dJ* (r,)p (C)R^R V M dr, d( ds (4.22c) 

- / / / Xse-^THriTHOiR^ + BlR^MdvdCds. (4.22d) 
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where 




(d v X)s + XsVsd v (f>. 



(4.23) 
(4.24) 



are symbols with a support of size ~ l/v^ an d enjoy uniform Coifman- Meyer bounds. Let us 
start with (|4.22ap : 



where we chose p > 6/5 sufficiently small; here we used Bernstein's inequality, as well as appealed 
to (|4.10p . The term (|4.22bj) can be treated similarly to (|4.22ap . thus we skip it. For (|4.22c|) . we 
estimate using Theorem 13. 131 that 



Here we have exploited the fact that Xs exhibits the improved Coifman-Meyer structure (in fact, 
on the frequencies < 1 where Xs is supported, the two notions of Coifman-Meyer symbols and 
improved Coifman-Meyer symbols coincide). We omit the details for the final term, (|4.22dD , 
because it is in fact a little simpler (no derivatives hit the profile). 

4.3.9. Weighted L 2 estimate for II ift < 1. Once again, it suffices to recapitulate the argument 
from the previous subsection, using the Sobolev embedding theorem where previously we had 
the dispersive estimates. 

4.3.10. H} estimates for J 7 " II. Multiplying by £ and evaluating the L? norm, we see that 



where m s is a Coifman-Meyer symbol supported on a ball of size l/\/s- 

For large time, the desired bound can be obtained with the help of the dispersive estimates, 
whereas we resort to Sobolev embedding for small time. We only illustrate the latter case: if 
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t < 1, 

< t ^\\e lsH f\\ L 4e isH f\\ L eds+ f ±\\ e isH f\\ L3 - S \\e isH f\\ L e ds 
» JO V s Jo V s 

< fJ-Jff ds <\\ff x . 
Jo V s " 

This completes the H} estimates for II, and the proof of Proposition 14.11 □ 
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